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Abstract. The qualitative behavior of a thermodynamically consistent two- 
phase Stefan problem with surface tension and with or without kinetic under- 
cooling is studied. It is shown that these problems generate local semifiows 
in well-defined state manifolds. If a solution does not exhibit singularities in 
a sense made precise below, it is proved that it exists globally in time and its 
orbit is relatively compact. In addition, stability and instability of equilibria 
is studied. In particular, it is shown that multiple spheres of the same radius 
are unstable, reminiscent of the onset of Ostwald ripening. 



1. Introduction 

Let f2 C K™ be a bounded domain of class C 2 , n > 2. £1 is occupied by a material 
that can undergo phase changes: at time t, phase i occupies the subdomain fli(t) of 
0, respectively, with i = 1,2. We assume that dili(t)ndfl — 0; this means that no 
boundary contact can occur. The closed compact hypersurface T(t) := dfli(t) C f2 
forms the interface between the phases. By the Stefan problem with surface tension 
we mean the following problem. 

Find a family of closed compact hypersurfaces {T(t)}t>o contained in Q and an 
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Here u(t) denotes the (absolute) temperature field, v(t) the outer normal field of 
Qi(i), V(t) the normal velocity of T(t), H(t) = H(T(t)) = -div r(t) i/(i)/(ra - 1) 
the mean curvature of T(t), and [t>] = «2|r(t) ~ w i|r(t) the jump of a quantity v 
across T(t). The sign of the mean curvature T-L is chosen to be negative at a point 
x € r if f2i PI B r (x) is convex for some sufficiently small r > 0. Thus if S7i is a ball 
of radius R then Ti = —1/R for its boundary T. 

Several quantities are derived from the free energies il>i(u) as follows: 

• ei(u) := ipi(u) + urji(u), the internal energy in phase i, 

• r)i(u) := — ?//(u), the entropy, 

• Ki(u) :— ti{u) = —uip'Ku) > 0, the heat capacity, 

• l(u) := = — u[?7(n)], the latent heat. 

Furthermore, di(u) > denotes the coefficient of heat conduction in Fourier's law, 
j(u) > the coefficient of kinetic undercooling, and a > the coefficient of surface 
tension. In the sequel we drop the index i, as there is no danger of confusion; we 
just keep in mind that the coefficients depend on the phases. 

The temperature is assumed to be continuous across the interface. However, 
the free energy and the conductivities depend on the respective phases, and hence 
the jumps [^(u)], [k(u)J, [^(tt)], [rf(w)| are in general non-zero at the interface. 
In this paper we assume that the coefficient of surface tension is constant. Then 
this model is consistent with the laws of thermodynamics. In fact, the total energy 
of the system is given by 

E(u,T)= [ e(u)dx + —^— / ads, (1.2) 
Jn\r n - 1 J T 

and by the transport and surface transport theorem we have for smooth solutions 
^E(«(t), r(t)) = - J^{\d(u)d v u] + \e{u)\V + aHV} ds 

= - j{{d{u)d v u\ - (l(u) - i{u)V))V} ds = 0, 

and thus, energy is conserved. Let us point out that it is essential that a > is 
constant, i.e. is independent of temperature. The reason for this lies in the fact 
that in case a = <j(u) depends on the temperature, the surface energy will be 
J r er(u) ds instead of J r ads, where er(u) — cr(u) + ur/r(u), rjr(u) = —a'(u), and 
one has to take into account the surface entropy J r r/r ds as well as balance of 
surface energy. This means that the Stefan law needs to be replaced by a dynamic 
boundary condition of the form 

n r {u)d t , n u - div r (dr(w)V r u) = \dd u u\ - (l{u) - j(u)V + l v {u)U)V, 

where dt, n denotes the time derivative in normal direction, nr{u) — ^'r( u ) an d 
lr{u) = ua'(u). We intend to study such more complex problems elsewhere and 
we restrict our attention to the case of constant a, here. 

The fifth equation in (|l.ip is usually called the Stefan law. It shows that energy is 
conserved across the interface. The fourth equation is the famous Gibbs-Thomson 
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law (with kinetic undercooling if j(u) > 0) which implies together with Stefan's 
law that entropy production on the interface is nonnegative if 7 > 0. In case 7 = 0, 
i.e. in the absence of kinetic undercooling, there is no entropy production on the 
interface, see below. Indeed, the total entropy of the system, given by 

$(u,T) = / r)(u)dx, (1.3) 
Jn\r 

satisfies 

^$(u(t),r(t))= [ \d{u)\Vu\ 2 dx- f -{{d{ U )d v ul+uHu)\V}ds 
at Jfi u Jr u 

= [ -^d(u)\Vu\ 2 dx + [ --f(u)V 2 ds > 0. 
Jn u z J T u 

In particular, the negative total entropy is a Lyapunov functional for problem 
(|l.ip . Even more, — $ is a strict Lyapunov functional in the sense that it is strictly 
decreasing along smooth solutions which are non-constant in time. Indeed, if at 
some time to > we have 

^$(«(t ),r(to)) = / \d(ti)\Vu\ 2 dx + f - 1 {u)V 2 ds = 0, 
dt J n u A J r u 

then Vit(io) = in Q and j(u(t ))V(t ) = on T(t )- This implies u(t a ) = const 
in Q, hence H(t ) = —\ijj{u{t n ))\/ a — const. Since ft is bounded, we may conclude 
that r(io) is a union of finitely many, say m, disjoint spheres of equal radius, i.e. 
(u(to), r(io)) is an equilibrium. Therefore, the limit sets of solutions in the state 
manifold SM. 1 , see (13.8l) - (|3.9p for a definition, are contained in the {ran + 1) 
dimensional manifold of equilibria 

£ = {(tt*, |J S R ,{xi)) :u* > 0, R* = cr/[V>(u*)], B R ,{xi) C ft, 

l<l<m 

andS Rt (x l )nS R ,(x k )=tD, I ? k}, (1.4) 

where S Rt (xi ) denotes the sphere with radius i?„ and center X\ . 

Another interesting observation is the following. Consider the critical points of 
the functional $(u,T) with constraint E(it, T) = E , say on C(Q) x MV. 2 (fl), 
see Section 3.1 for the definition of A4H 2 (n). Then by the method of Lagrange 
multipliers, there is \x € K such that at a critical point (ii*,L») we have 

$'(u*,L») + mE'(w*,L„) =0. (1.5) 

The derivatives of the functionals are given by 

(&(u,r)\(v,h)) = W(u)\v) L2m - (Mu)}\h) L2{r) , 

and 

(E'(u,T)\(v,h)) = (e'(u)\v) L2(n) - (le(u)]+aH(T)\h) L2{ ry 
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Setting first h = and varying v in (II. 5j) wc obtain rf (it*) + /ie'(it*) = in fi, and 
then varying h we get 

fo(u,)] + ^([e(u,)l + <TW(r,)) = on r„. 

The relations r?(u) = —ip'(u) and e(u) = i/)(u)— uip'(u) imply = — ^>"(u*)(l+/X'ii*), 
and this shows that w* = —1/fJ, is constant in fi, since k(u) = —uip"(u) > for all 
u > by assumption. This further implies [?/>(«*)] + a'HiT*) — 0, i.e. the Gibbs- 
Thomson relation. Since it* is constant we see that HiT*) is constant, hence T* is 
a sphere whenever connected, and a union of finitely many disjoint spheres of equal 
size otherwise. Thus the critical points of the entropy functional for prescribed 
energy are precisely the equilibria of the problem. 

Going further, suppose we have an equilibrium e* := (it*,r„) where the total 
entropy has a local maximum, w.r.t. the constraint E = Eo constant. Then D := 
[$ + /iE]"(e*) is negative semi-definite on the kernel of E'(e*), where fx is the fixed 
Lagrange multiplier found above. The kernel of E'(e) is given by the identity 

k(u)v dx — 

which at equilibrium yields 

n*v dx + / l*hds = 0, (1-6) 
n Jr, 

where k* := rt(ti*) and l* :— l(u*). On the other hand, a straightforward calcula- 
tion yields with z — (v,h) 



— j k»v 2 dx — a f h-7i'(T*)hds 



(1.7) 



As k* and a are positive, we see that the form (T>z\z) is negative semi-definite as 
soon as ^'(r*) is negative semi-definite. We have 

H'(r*) = l/i?2 + (l/(n-l))Ar,, 

where Ap, denotes the Laplace-Beltrami operator on r* and i?* means the radius 
of the equilibrium sphere. To derive necessary conditions for an equilibrium e* to 
be a local maximum of entropy, we consider two cases. 

1. Suppose that T* is not connected, i.e. r* is a union of m disjoint spheres Tj. 
Set v = 0, and let h = h/. ^ be constant on r* (and zero on the remaining 
components of T*) with ^2 k hf. = 0. Then the constraint (|1.6[) holds, and 

(Vz\z) = (a/u*Rl)(\T*\/m)^2h 2 k > 0, 

k 

hence T> cannot be negative semi-definite in this case. Thus if e* is an equilibrium 
with maximal total entropy, then r* must be connected, and hence both phases 
are connected. 
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2. Assume that T* is connected. Setting v = Z*/(«*|l)n and h = — l/|r*|, we see 
that T> negative semi-definite on the kernel of E'(e*) implies the condition 

c * - * (L8) 

This is exactly the stability condition found in Theorem 14.51 
In summary, we have shown that: 

• The equilibria of (11.11) are precisely the critical points of the entropy func- 
tional with prescribed energy. 

• The entropy functional with prescribed energy does not have local maxima 
in case I\ is not connected. 

• In case T* is connected, a necessary condition for a critical point (tt*,r*) 
to be a local maximum of the entropy functional with prescribed energy 
is inequality (ll.8[) . 

It will be shown in Theorems 14.51 and 15.21 below that 

• («*, r») <E £ is stable if r» is connected and £* < 1< 

• The latter is exactly the case if the reduced energy functional, 

[u i-> (f(u) = E(u, S R{u) (x ))}, R(u) = a/{^{u)\, 

has a strictly negative derivative at u*. 

• Any solution starting in a neighborhood of a stable equilibrium exists 
globally and converges to another stable equilibrium exponentially fast. 

• (it*, r„) g £ is always unstable if is disconnected, or if (* > 1. 

Hence multiple spheres (of the same radius) are always unstable for (|1.1[) . This 
situation is reminiscent of the onset of Ostwald ripening. 

The authors in pQ also consider the case with the Gibbs-Thomson law replaced by 

li>(u)] + aU = 7(tt)y-divr[a(it)Vr(V7ii)] -udivr[/3(u)VrV], 

where a, /3 > may also depend on the temperature. The modified Stefan law in 
this case reads 

{d{u)d v u\ = (l(u) -j(u)V + div r [a(u)V r (V r /u)] + udiv r [/3(u)V r V r ]) K 
The resulting entropy production then becomes 

iU(u(t),r(t))= / \d{u)\Vu\ 2 dx+ f - 7 (u)V 2 ds 
dt Jq u Jt u 



f J [a(u)\\7 r (V/u)\ 2 + (3(u)\V r V\ 2 ]\ds > 0. 



Although interesting from an analytic point of view, these models will not be 
considered here. 
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Now we want to relate problem (jl.lj) to the Stefan problems that have been studied 
in the mathematical literature so far. For this purpose we linearize h(u) := [y>(it)] 
near the melting temperature u m defined by h(u m ) = 0. Then for the relative 
temperature v = u — u m we have h(u) ~ h'(u m )v, hence with l m = l(u m ) and 
7m = 7(1*171)7 the Gibbs-Thomson law approximately becomes 

{lm/u m )v + (TV. = 7 m V. (1.9) 

This is the classical Gibbs-Thomson law with kinetic undercooling. Similarly, as- 
suming that u is close to u m and V is small, the Stefan law becomes approximately 

\dd v u\ = l m V. (1.10) 

In the literature it is frequently assumed that the heat capacities are constant and 
equal in the phases. Then we necessarily have 

0=14 = [e»] = -t*[V»"(«)l> 

which implies that the function h(u) — [ip{u)l is linear, i.e. h(u) = ho + h\u, and 
then l(u) = h\u. The melting temperature is here given by < u m — —ho/h\. 

If the heat capacities Ki are constant in the phases but not necessarily equal, 
the internal energies depend linearly on the temperature and the free energies are 
of the form ipi(u) = a, + biU — Kjitlnti, hence h(u) = a + /3u — Suln(u), with 
constants a, (3,6 € M. Concerning existence of equilibria, these special cases will 
be discussed in more detail in Section 4. 

Let us now briefly discuss the literature concerning problem (jl.ip and related 
problems. We refer to pQ and [22] for information on the modeling aspects of (jl.lll . 
see also Gurtin [T8l H"9] for earlier work. Additional information for the classical 
Gibbs-Thomson law (H) can be found in [H [H [HI [HI EH EH [38] . 

Existence and regularity results for the Stefan problem with the classical Gibbs- 
Thomson law u = aT-L and the classical Stefan law (|1.10[) in case K\ = K2 can be 
found in [H EH Hi Hi H3 [351 see [H] for a detailed discussion. The Stefan 
problem with the Gibbs-Thomson law (|1.9p and the classical Stefan law (jl.lOp in 
case Ki = k 2 has been studied in [71 133 [351 [37], see also [53] for the one-phase 
case. 

If «i = K2 = then we obtain a thermodynamically consistent quasi-stationary 
approximation of the Stefan problem with surface tension (and kinetic undercool- 
ing). Existence, uniqueness, regularity, and global existence of solutions for the 
quasi-stationary approximation with the Gibbs-Thomson law u — aH and the clas- 
sical Stefan law (Ti~TU)) has been studied in [U [H O [13 [HI HU H3J [H] • Existence 
and global existence of classical solutions for the quasi-stationary approximation 
with (|1.9p and the classical Stefan law (|1.10[) has been investigated in [39[ [23] . 

We refer to [3U] for a detailed discussion of the literature of the classical Stefan 
problem with K\ = Ki, o = 7 = and (| 1 . 1 0|) . 

It appears that this manuscript is the first work to provide analytical results 
for the thermodynamically consistent Stefan problem (II. 
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The plan for this paper is as follows. In Section 2 we first transform the problem 
to a domain with fixed interface, employing a Hanzawa transform |20j . This is the 
so-called direct mapping approach. Section 3 is devoted to results on local well- 
posedness for problem , based on the approach in [15] and [9] . We introduce 
the state manifolds SA4 1 for (jl.ll) and we consider the local semiflow in SAi^ 
generated by problem . In Section 4 we discuss equilibria and their linear sta- 
bility properties. Here we rely on previous work of the authors [32] ■ In Section 5 
we establish the corresponding stability properties for the nonlinear problem, em- 
ploying the generalized principle of linearized stability, extending the results of |33] 
to the situation considered here. The main result of this section shows convergence 
of solutions to an equilibrium which start out near stable equilibria. Moreover, we 
give a rigorous proof of the instability result. Of ultimate importance is the Lya- 
punov functional for which is given by the negative total entropy — <E>(zi,r). 

It takes bounded global-in-time solutions to the set of equilibria, and then by the 
results of Section 5 and compactness of the orbits, such solutions must converge 
towards an equilibrium in the topology of the state manifold SM.^ provided they 
come close to a stable equilibrium. 

Our analysis is carried out in the framework of L p -spaces, with n + 2 < p < oo. 
We expect that it would be enough to require (n + 2)/2 < p < oo (so unfortunately 
p > 2 even in 2DI), but for the sake of simplicity we restrict ourselves here to the 
stronger assumption p > n + 2. We also expect that a similar analysis can be 
obtained in the framework of the small Holder spaces h a , which would, though, 
require higher order compatibility conditions. 

2. Transformation to a Fixed Interface 

Let ft C R™ be a bounded domain with boundary 9ft of class C 2 , and suppose 
r C ft is a closed hypersurface of class C 2 , i.e. a C 2 -manifold which is the boundary 
of a bounded domain fti C ft. We then set ft2 = ft \ fti. Note that while ft2 
is connected, fti may be disconnected. However, fti consists of finitely many 
components only, as <9fti = T by assumption is a manifold, at least of class C 2 . 
Recall that the second order bundle of T is given by 

^r:={(p,i<r(p),Lr(p)) : P G T}. 

Note that the Weingarten map Lr (also called the shape operator, or the second 
fundamental tensor) is defined by 

£r(p) = -VrtrO), per, 

where Vr denotes the surface gradient on T. The eigenvalues Kj{p) of Lr(p) are 
the principal curvatures of T at p € l\ and we have |Lr(f>)| = max, |rej(p)|. The 
mean curvature 'Hr(p) is given by 

n-i 

(n - l)H r (p) = ^2 Kj(p) = trirO) = -div r ^r(p), 
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where divr means surface divergence. Recall also that the Hausdorff distance dn 
between the two closed subsets A,B C M. m is defined by 

da (A, B) := max { supdist(a, B), sup dist(6, A)}. 

aeA beB 

Then we may approximate T by a real analytic hypersurface £ (or merely £ € C 3 ), 
in the sense that the Hausdorff distance of the second order bundles of T and £ 
is as small as we want. More precisely, for each 77 > there is a real analytic 
closed hypersurface such that df/(A/" 2 £, J\f 2 T) < t). If 77 > is small enough, then 
E bounds a domain fif with flf C O, and we set fif = SI \ Clf . 

It is well known that such a hypersurface £ admits a tubular neighborhood, 
which means that there is a > such that the map 

A : £ x (-a, a) -> M" 
A(p, r) := p + rvT.{p) 

is a diffeomorphism from £ x (—a, a) onto 7\L(A). The inverse 

A" 1 : TZ{A) 1— > £ x (-a, a) 

of this map is conveniently decomposed as 

A _1 (x) = (n(x),d E (x)), jceft(A). 

Here n(x) means the nonlinear orthogonal projection of x to E and <fe(x) the 
signed distance from x to E; so |ds(x)| = dist(x, E) and ds(x) < iff x G Of. In 
particular we have 1Z(A) = {x e R" : dist(x, E) < a}. 

On the one hand, a is determined by the curvatures of E, i.e. we must have 

< a < min {l/|/tj(p)| :j = l,...,n-l,peE}, 

where mean the principal curvatures of E at p € E. But on the other hand, 

a is also connected to the topology of E, which can be expressed as follows. Since 
E is a compact (smooth) manifold of dimension n — 1 it satisfies a (interior and 
exterior) ball condition, which means that there is a radius > such that for 
each point p e E there are Xj £ itj, j — 1,2, such that B rs (xj) C Slf, and 
B rs (xj) n E = {p}. Choosing maximal, we then must also have a < r%. In the 
sequel we fix 

a = iminirs, - — j = 1, . . . , n - 1, p e E 1 . 
For later use we note that the derivatives of n(x) and (is(x) are given by 

vd E (x) = i^(n(aO), n'(x) = M (d s (x),n(x))p s (n(x)), 

where Ps(p) = I — vy<{p) ® v t.{p) denotes the orthogonal projection onto the 
tangent space T p E of E at p e E, and M (r,p) = (I — rL^(p))^ 1 . Note that 
\M (r,p)\ < 1/(1 - r|L E (p)|) < 2 for all |r| < a and p e E. 
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Setting r = r(i), we may use the map A to parameterize the unknown free 
boundary T(t) over E by means of a height function p(t,p) via 

T(t):]p^p + p{t t p)vE(p)], peE,t>0, 

for small t > 0, at least. Extend this diffcomorphism to all of Ct by means of 

E p (t, x)=x + x(ds(x)/a)p(t, IL(x))vx (1%)) =: x + 6 p (t, x). 

Here x denotes a suitable cut-off function; more precisely, \ £ £>(K0, < x < lj 
X(r) = 1 for |r| < 1/3, and x(r) = for |r| > 2/3. Note that E p (t,x) = a; for 
|d s (ar)| > 2a/3, and 

H-^t.a;) =ar-p(t,n(:r)K(n(ar)) for \d s (x)\ < a/3. 

In particular, 

S~ 1 (t, x) = a; — p(t, x)uy,{x) for ieE. 

Setting v(t, x) — u(t, E p (t, x)), or u(t,x) — v(t, E^ 1 (t, x)) we have this way trans- 
formed the time varying regions fi \ T(t) to the fixed domain Q \ E. This is the 
direct mapping method, also called Hanzawa transformation. 

By means of this transformation, we obtain the following transformed problem. 

n(v)dtv + A(v, p)v = n(v)lZ{p)v 
d Vn v = 

H = o 

W(v)] + aH(p) = j(v)(3(p)d tP 
{l(v) - 1 (v)P(p)d tP }P(p)dtP + B(v, p)v = 
v(0) = v , p(Q) = p Q . 

Here A(v,p) and B(v,p) denote the transformations of — div(dV) and — [cZ9„], 
respectively. Moreover, H(p) means the mean curvature of T, j3(p) — {v?\vr(p)), 
the term f3(p)dtp represents the normal velocity V, and 

lZ(p)v = d t v - d t u o E p . 

The system (|2.ip is a quasi-linear parabolic problem on the domain ft with fixed 
interface S C O with a dynamic boundary condition, namely the fifth equation 
which describes the evolution of the interface T(t). 

To elaborate on the structure of this problem in more detail, we calculate 

% = I + 8> p , E'- 1 =/-[/ + & p \~\ =: / - M 1 (p) T 

and 

Vu°E p = [(E- 1 )' 7 o E p ]Vv = (I - Mi(p))V«, 
and for a vector field q — q o E p 

(V\q) o E p = ([(E^ f o E p ]V\q) = ((I - Mx(p))V|g). 



m 

on 
on 
on 
on 



dQ 
E 
E 
E 



(2.1) 



10 



J. PRUSS, G. SIMONETT, AND R. ZACHER 



Further wc have 



d t u o Z p = d t v - (Vu o E p \d t E p ) = d t v - ([(E p 1 )' T o E p ]Vi;|0 t S p ) 
= ft«-(Vt;|[/ + ^]- 1 ^) ) 

hence 

7e(p> = (V«|[/ + ^]- 1 9 t ^). 
With the Weingarten map L E = — Vji/j we have 

= P{p){vn - a(p)), a(p) = M (p)V E p, 
M (p) = (/ - pL^)-\ p{p) = (1 + la(p)! 2 )- 1 / 2 , 

and 

V = (d t E\ur) = {vY,\v T {p))d t p = (3(p)d t p. 
Employing this notation leads to 9' p = for |d E (x)| > 2a/ 3 and 

6' p (t,x) = ^x'(d^(x)/a)p(t,U(x))^(U(x))®^(U(x)) 

+ X (d E (a:)/a)[i^(n(a;))(8i M (d E (a;))VEP(t,n(a;))] 
-x(ds(a:)/a)p(t,n(a;))LE(n(a;))Mo(d E (!c))ft(n(x)) 

for < |<fe(x)| < 2a/3. In particular, for x e £ we have 

p (t,x) = ^ E (x) <g> V S /3(i,i) - p(t, x)L E (x)P E (x), 

and 

(6> p ) T (t, x) = V E p(i, a;) <8> ws(a;) - p(t, x)L s (x), 

since Le(i) is symmetric and has range in T^E. Therefore, [/ + Q' p \ is boundedly 
invertible, if p and V E/ o are sufficiently small, and 

For the mean curvature "H(p) we have 

(n - l)H(p) = /3(p){tr[M (p)(L s + V E a(p))] - /3 2 (p)(M (p)a(p)| [V E a(p)]a(p))} ; 

an expression involving second order derivatives of p only linearly. Its linearization 
at p — is given by 

(n- 1)W(0) =trL| + A E . 

Here A E denotes the Laplace-Beltrami operator on E. The operator B(v,p) be- 
comes 

B(v,p)v = -ld{u)d v u\ o E p = -([d(v)(I - Mi(p))V«]|i/r) 
= ([<*(<;) (J - M x (p))Vv\\w - a( P )) 

= -f3(p)ld(v)d vs v] + /3(j>)(ld(v)Vv]\(I-M 1 (p)) T a(p)), 
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since Mj(p)vj: — 0, and finally 

A(v, p)v = - div(d(u)Vw) o E p = -((/ - Mi(p))V|d(u)(7 - M^Vv) 
= - d(v)Av + d(v)[M!(p) + Mj{p) - Mi{p)Mj{p)] : V 2 v 
- d'(v)\(I - M 1 (p))Vw| 2 + d(v)((I - Mi(p)) : VMi(p)|V«). 

We recall that for matrices A, B E R nXTi , A : B = Yn,j=i <^ijh = tr (AB T ) 
denotes the inner product. 

Obviously, the leading part of A(v,p)v is — d(v)Av, while the leading part of 
B(v,p)v is — [3{p)\d(v)d Ul ,v\, as Mi(0) = and a(0) — 0; recall that we may as- 
sume p small in the C 2 -norm. It is important to recognize the quasilinear structure 
of (j2.1[) : derivatives of highest order only appear linearly in each of the equations. 

3. Local Well-Posedness 

The basic result for local well-posedness in the absence of kinetic undercooling 
in an L^-setting is the following. 

Theorem 3.1. (7 = 0). Let p > n + 2, 7 = 0, a > 0. Suppose tp € C 3 (0,oo), 
d G C 2 (0,oo) such that 

k(u) = -uip"(u) > 0, d(u) > 0, ue{0, 00). 

Assume the regularity conditions 

u ew*- 2 ^(n\r )nc({i), u >o, r ew^ 3 /p, 

the compatibility conditions 

[V>M] + aU{T Q ) = 0, {d{u a )d VT< u \ € W^- 6 /*> (r ), 
and the well-posedness condition 

K u o) 7^ on Tq. 

Then there exists a unique L p - solution for the Stefan problem with surface tension 
on some possibly small but nontrivial time interval J = [0, r]. 

Here the notation To € Wp~ 3 ^ p means that To is a C 2 -manifold, such that its 
(outer) normal field vy is of class Wp 3 ^ p (ro). Therefore the Weingarten tensor 

Lr a = — Vr iT of r belongs to W%~ 3 ^ p (Tq) which embeds into C 1+Q (r ), with 
a = 1 — (n + 2)/p > since p > n + 2 by assumption. For the same reason, we 
also have uq g C 1+Q (f2), and Vb € C 2a (ro). The notion L p -solution means that 
(w,r) is obtained as the push-forward of an L p -solution (v,p) of the transformed 
problem (|2.1I) . This class will be discussed below. 

There is an analogous result in the presence of kinetic undercooling which reads 
as follows. 
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Theorem 3.2. (7 > 0). Let p > n + 2, a > 0, and suppose ifj,j & C 3 (0, 00), 
d G C 2 (0, 00) such that 

K (u) = -wip"(u) > 0, d(u) > 0, 7(u) > 0, 116 (0,oo). 

Assume the regularity conditions 

u a e w 2 ~^ p (n \ v ) n C(O), u > o, r e , 

and the compatibility condition 

([VM] + o-H(T ))(l(u ) - [V»(«o)] - o^(T )) = 7 MNMftuo]- 

TTien i/iere exists a unique L p -solution of the Stefan problem with surface tension 
and kinetic undercooling (11. lj) on some possibly small but nontrivial time interval 
J=[0,t}. 

Proof of Theorems GO] and EH} 

(i) Direct mapping method: Hanzawa transformation. 

As explained in the previous section, we employ a Hanzawa transformation and 
study the resulting problem (|2.1[) on the domain 51 with fixed interface E. 

In case 7 = 0, for the L p -theory, the solution of the transformed problem will 
belong to the class 

v G Hp(J; L p (Sl)) n L p (J; H 2 (n \ E)) C(J; W^^n \ E)), 
p G ^-^(J; # 2 (E)) n L p (J; % 4 -Vf (S)) ^ C(J; (E)), 

ftp G W p 1 /2-l/2p ( J. Lj>(E)) n Lp( J. W 2-2/ P(s)) ^ C(J . w 2-6/„ (E))) (31) 

see [15] for a proof of the last two embeddings in the case E = W 1 . 
If 7 > we have moreover 

p G W^ 2 P(J; L P (E)) n L P (J; W^{T)) ^ C 1 (J; W 2 p -^(Z)). 

Note that in both cases, v G C(J x O), G C(J; C^fij)), j = 1,2. Moreover, 
pe C(J;C 3 (E)) and 

ftp G C*(J; C(E)) in case 7 = 0, ftp G C(J; C^E)) in case 7 > 0. 

We set 

Ex (J) := {v G H^J-L p (nj) n L p (J;H 2 (n \ E) : [«] = 0, ft„u = 0}, 

E 2 (J) := Wp^ 2 ~ 1 S 2p (J; L P (E)) n W^-V^J; H p 2 (E)) n L P (J; W^~Vf (E)), 7 = 0, 

Ea( J) := % 2 " 1/2p ( J; L P (E)) n L p { J; <~ 1/p (E)), 7 > 0, 

E(J) :=Ei(J) x E 2 (J), 
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i.e. E( J) denotes the solution space. Similarly, we define 
F X (J) :=L P (J;L P (Q)), 

F 2 (J) := W^ 1 - 1 / 2 f(J;i p (I]))nL p (J;W^- 1 ^(E)), 

F 3 (J) := Wpt 2 ~ 1 t 2p {J; L p (E)) n L p (J; (£)), 

F(J) :=F 1 (J)xF 2 (J)xF 3 (J), 

i.e. F( J) means the space of data. A left subscript zero means vanishing time trace 
at t = 0, whenever it exists. So for example 

E 2 (J) = {pe E 2 (J) : p(0) - 9 4 p(0) = 0} 

whenever p > 3. Employing the calculations in Section 2 and splitting into the 
principal linear part and a nonlinear part, we arrive at the following formulation 
of problem (|2~Tj) . 

no(x)dtv — do(x)Av 
d vsi v 

H 

x)v + cr A^p - 7i (t, x)d t p 
l (x)d t p - {do{x)d v v\ 
u(0) = v , p(0) 

Here 

k (x) = k(v (x)), d (x) = d(v (x)), l (x) = l(v (x)), a Q = 
h(t, •) = [V>'(e As W)], 7i (*, ') = 7(e As W), 

where Uos means the restriction of Vq to S. Note that Ko,<io € Wp~ 2 ^ p (tt \ £), 
hence these functions are in C 1 (f2j), j = 1,2. Recall that <i and n may be different 
in different phases. Further we have lo € Wp 2 3 ^(S) which implies Iq € C 1 (E). 
This is good enough for the space F 3 ( J), as C 1 -functions are pointwise multipliers 
for F 3 (J), but it is not good enough for F 2 (J). For this reason, we need to define 
the extension Vb ■— e As *wos- This function as well as l\ and 71 belong to F 2 (J), 
hence are pointwise multipliers for this space, as F 2 (J) and ¥s(J) are Banach 
algebras w.r.t. pointwise multiplication, as p > n + 2. 
The nonlinearities F, G, and H are defined as follows. 

F(v, p) = (k - n{v))d t v + (d(v) - d Q )Av + d{v)M 2 {p) : V 2 v 

- d'(v)\(I - Ahip^Vv] 2 + d{v)(M 3 (p)\Vv) + k(v)K( P )v, 
G(v, p) = + aU{p)) + hv + aoA^p + ( 7 (»)/3(p) - ji)d t p, (3.3) 

H(v,p) = l(d(v) - d )d v v] + (l - l(v))d tP + ([d(v)Vv]\M 4 (p)Vxp) 
+ 1 (v)(3(p)(d t p) 2 . 



F(v,p) in Q\S 

on dfl 

0, on £ 

G(v,p) on £ 

H (v, p) on £ 
Pa- 



(3.2) 
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Here we have set 

M 2 (p) = M x (p) + AlJ(p) - M 1 (p)M 1 T ( /) ), 
M 3 (p) = (/-M 1 (/>)) :VMi(p), 
M 4 (p) = (/-M 1 (p)) T M (p). 

(ii) Maximal regularity of the principal linearized problem. 

First we consider the linear problem defined by the left hand side of (|3.2p . 



K Q (x)d t v 


— do(x) At; 


= / 


in 


n\s 






= 


on 






M 


= 


on 


E 


h(t,x)v + f7 A s p - 


j 1 (t 1 x)d t p 


= .9 


on 


E 


l {x)d t p- 


[d (x)d v vj 


= h 


on 


E 


«(0) 


= v , p(0) 


= Po- 







(3.4) 



This inhomogeneous problem can be solved with maximal regularity; see Escher, 
Priiss, and Simonett [15] for the constant coefficient half-space case with 7 = 0, 
and Denk, Priiss, and Zacher [3] for the general one-phase case. 

Theorem 3.3. (7 = 0). Let p > n + 2. er > 0, 7 = 0. Suppose kq G C(flj) and 
d E C^Clj), j = 1,2, K ,d > on Q, l G Wp 2 ~ 6/p (E), and let 

h G W p 1 - 1 /2 P (J;i p (E))ni p (J;W p 2 - 1 /f(E)) 

such that IqI\ > on J x E, where J — [0,io] is a finite time interval. Then there 
is a unique solution z := (v, p) G E(J) 0/ i/ and on/y if the data (/,<?, /1) and 

z :=(v ,po) satisfy 

(f,g,h)€¥(J), ^[^(^ElnCix^lE), 

and the compatibility conditions 

h(0)v + CTo A s po - .9(0), MO) + [<hd v vo] G W*- 6/p (Z). 

The solution map [(/, g, h, Zq) > z = (v,p)] is continuous between the correspond- 
ing spaces. 

Proof. In the one-phase case this result is proved in 9, Example 3.4] . Therefore, we 
only indicate the necessary modifications for the two-phase case. The localization 
procedure can be carried out in the same way as in the one-phase case [9], hence 
we only need to consider the following model problem with constant coefficients 
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where the interface is flat: 

KodtV — dAv = f in 

[vj = on 

< l\v + (ToAp = g on 

lodtp - \dd v v\ =h on 

t v(0) = v , p(0) = p . 

Here R™ = l"" 1 x (R\ {0}), and R™ -1 is identified with R n_1 x {0}. Reflecting 
the lower half-plane to the upper, this becomes a problem of the form studied 
in [9]. As in Example 3.4 of that paper it is not difficult to verify the necessary 
Lopatinskii-Shapiro conditions. Then Theorems 2.1 and 2.2 of [5] can be applied, 
proving the assertion for the model problem. □ 

Remark 3.4. One might wonder where the somewhat unexpected compatibility 
condition h(0) + [do^'Wo] £ Wp~ 6 ^ '(E) in the case 7 = comes from. To allude 
this, note that 

(h(0) + [dod v vo])/lo = ftp(0) 
is the trace of dtp at time t = 0. But by the embedding fl3.1[) this implies that 
(h(0) + [doduVoD/lo £ Wp 2_6/p (S), which in turn enforces h(0) + [do<%w ] £ 
W p 2 ~ 6/p (E). 

The main result for problem (|3.4I) for 7 > is the following theorem. 

Theorem 3.5. (7 > 0). Let p > n + 2, a > 0. Suppose kq £ C(Ctj) and 
d a £ C 1 ^), j = 1,2, n a ,d a > on fl, l £ C^E), and Zei 

71, Zi € iy p 1 - 1 / 2p (J;L p (S))nL p (J;^ p 2 - 1 / p (E)), 

smc/i that 71 > on J x S, where J = [0,to] is a finite time interval. Then there 
is a unique solution z :— (v, p) £ E( J) 0/ \3.4\ ) if and only if the data (/, g, h) and 
z := («o,Po) 

(f,g,h)£F(J), zoe^t^sjnc^x^fs), 

and i/ie compatibility conditions 

(loh{0)v \ s + ZofToAsPo - 7i(0)[dd„fo] = 71 (0)fc(0) + «off(0). 

TTie solution map [{f,g,h,Zo) z — (v,p)] is continuous between the correspond- 
ing spaces. 

Proof. The proof of this result is much simpler than for the case 7 = 0. We could 
follow the strategy in the proof of Theorem l3.31 employing the methods in [5] once 
more. However, here we want to give a more direct argument that uses the fact 
that the term l$dtp is of lower order in case 71 > 0. For this purpose, suppose 
v y. '■= v\ % is known. Consider the problem 

jidtp - o- A s p = Ixv-z - g, t£j, p(0) = p a . 



R 
R'' 
R' 



1(5 



J. PRUSS, G. SIMONETT, AND R. ZACHER 



Since the Laplace-Beltrami operator is strongly elliptic, we can solve this problem 
with maximal regularity to obtain p in the proper regularity class. Then we solve 
the transmission problem 



' KodtV — dgAv 


= f 


in 


n\s 


dvnV 


= 


on 


dCl 


M 


= 


on 


E 


-\dod v v\ 


= h — l Q d t p 


on 


E 


v(0) 


= va- 







Finally, we take the trace of v to obtain an equation for vs of the form 

Us = Tv% + w, 

where w is determined by the data alone, and T is a compact operator from F2 
into itself. Here compactness follows from the compact embedding F2 F3, i.e. 
from the regularity of dtp which is higher than needed to solve the transmission 
problem. Thus / — T is a Fredholm operator with index zero, hence invertiblc 
since it is injective by causality. This proves the sufficiency of the conditions on 
the data. Necessity is a consequence of trace theory. □ 

Remark 3.6. It is interesting to take a look at the boundary symbol of the linear 
problem. It is of the form 

«(A,0 = A; 2 + (A 7 + <7o|£| 2 )h/A«i + dilCI 2 + + d 2 \tf}. 

Here A £ C+ denotes the covariable of time t, and £ £ R™ -1 that of the tangential 
space variable x' £ M 71 . This symbol is invertible for large A, provided 7 > or 
I 7^ 0. Note that in case 7 > this is a parabolic symbol of order 3/2 in time t 
and of order 3 in the space variables x. The term XI 2 is of lower order, thus I does 
not affect well-posedness. On the other hand, for I = and 7 = the boundary 
symbol is ill-posed, since it admits the zeros (A, 0) with arbitrarily large Re A. If 
7 = and I ^ 0, then it is well-posed. Note that in this case we have order I in 
time, 3 in space, but also the mixed regularity 1/2 in time and 2 in space. 

(hi) Reduction to zero initial values. 

It is convenient to reduce the problem to zero initial data and inhomogeneities 
with vanishing time trace. This can be achieved as follows. We solve the linear 
problem (|3.4|) with initial data vo, po and inhomogenities 

/ = 0, g{t) = e Ast G{vo, p ), h{t) = e Ast pi with p x = H(v Q , p ). 

Since the Laplace-Beltrami operator As has maximal L p -regularity, the fact that 
G(v ,po) £ W / p 2_3/p (E) implies g £ F 2 (J). Similarly, /igF 3 since p x £ Wrl" 3/p (E). 
The compatibility conditions yield [do^wo] + Pi <= 

Wp~ 6/p (E). Therefore, the 
linear problem has a unique solution z* := with maximal regularity £ 
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(3.5) 



E(J). Then we set v — v ~ v„, p = p — p*, and obtain the following problem for 
z = {v,p). 

Ko(x)d t v — do(x)Av = F(v + v*, p + p*) in fi\E 

d m v = on dVL 

\v\ = on E 

h(t,x)v + a A^p - 7i(t,x)d t p = G(v, p;v*, p*) on E 

l Q (x)d t p - ldo(x)d v v} = £f(w,p;u*,p*) on E 
5(0) = 0, p(0) = 0. 

Here we have set 

G(v,p;v*,p*) = G(v + v*,p + p*) - e Ast G(v Q ,p ), 
H(v, p; i>*, p*) = #(u + u», p + p*) - e Ast H(v 0l po). 

Note that G(0, 0; Uoj Po) = H(0, 0; uo, Po) = by construction, which ensures time 
trace zero at t = 0. 



(iv) Solution of the nonlinear problem 

We first concentrate on the case 7 = 0, and rewrite problem (I3.5[) in abstract form 

as 

where L : E(0,t ) -> oF(0,t o ), defined by 

hz = (n d t v - d Av,hv + a A^p, l d t p - [dod v vf), 
is an isomorphism by Theorem 13.31 The nonlinearity 

N: o E(0,t o ) xE(0,toWoF(0,i o ), 

given by the right hand side of ([33]), is of class C 1 , since the coefficient functions 
satisfy k £ C 1 , d, I 6 C 2 , tp & C 3 , and by virtue of the embeddings 

Ei (J) <^> C(J x n) n CfJjC 1 ^))) E 2(J) ^ C*(J; C* 3 (E)) n C7 1 (J;C*(E)). 

Observe that the constants in these embeddings blow up as to — ¥ 0, however, they 
are uniform in to if one considers the space oE(J)! 

We want to apply the contraction mapping principle. For this purpose we 
consider a closed ball B#(0) C oE(0, r), where the radius R > and the final time 
t € (0, fo] are & t our disposal. We rewrite the abstract equation hz = N(z, z*) as 
the fixed point equation 

z = L -1 N(z, z.) =:T(z), zeM R (0). 

Since we are working in an L p -setting, by choosing r = t(R) small enough we can 
assure that 

||T(0)|| E(0 , T) = [L-^O,**)!^) < R/2. 
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On the other hand, we have 

|T(zi) - T(z 2 )| E (o,r) < |L _1 |i 3 ( F(0,r),oE(0,r)) x 

X Sup IN / (f,Z*)|| B ( E(0,T), F(0 1 T))I«l - «2lE(0,T)) 

II z II e(o,t) <-R 

hence T(Bfl(0)) C 1Br(0) and T is a strict contraction, provided we have 

||b( d F(0,t), o E(0,t)) Sup ||N'(z, Z*) |b( o E(0,t), o F(0,t)) < 1/2. 

ll z ll E(o, T )<« 

For this we observe that 

||b( o F(0,t), o E(0,t)) < le( o F(0,to),oE(0,t o )) = : C M < OO 

is uniform in t G (0, to), since we have vanishing time traces at t = 0. So it remains 
to estimate the Frechet-derivative of N on the ball B_r(0) C oE(0,t). This is the 
content of the next proposition, which also covers the case 7 > 0. 

Proposition 3.7. Let p > n + 2, a E M, and suppose ^,7 € C 3 (0, 00) and 
d E C 2 (0,oo). Then N : E(0,*o) X E(0,*o) -> oF(0, t ) *s continuously Frechet- 
differentiable. There is r\ > swc/i that for a given z* £ E(0,to) with |polc 2 (s) < V> 
there are continuous functions a(R) > and /3(r) > with a(0) = /?(0) = 0, such 
that 

||N'(z + z*)||b(qE(o,t), f(o,t)) < a(ii) + /3(r), zeB H c E(0, r). 

Proof. We may proceed similarly as in [151 Section 7] , where the interface is a graph 
over R n_1 . The additional terms which arise by considering a general geometry 
are either of lower order or of the form M(v, p)Vy,p where M(v,p) is of highest 
order (see (13.31) ). but can be controlled by ensuring that is sufficiently small. 
The additional terms due to the presence of 7 are of highest order, but small. □ 

So choosing first R > and then r > small enough, T will be a self-map and a 
strict contraction on 18^(0). Concluding, the contraction mapping principle yields 
a unique fixed point z = z(z^) £ 1Br(0) C oE(0,t), hence z — z* + z(z») is the 
unique solution of (|3.2j) . i.e. of (|2.ip . 

The proof in case 7 > is similar, employing now Theorem 13.51 

Remark 3.8. The assumption p > n + 2 simplifies many arguments since F 2 (J) 
as well as J) are Banach algebras and Vu E BC(J x 57). If we merely assume 
p > (n + 2)/2 then F 2 (J) is still a Banach algebra, but Fa(J) is not, and Vv may 
not be bounded anymore. This leads to much more involved estimates for the 
nonlinearities. 

Local semiflows. We denote by MH 2 ^) the closed C 2 -hypersurfaces contained 
in 51. It can be shown that A4H 2 (fl) is a C 2 -manifold: the charts are the param- 
ctcrizations over a given hypersurface E according to Section 2, and the tangent 
space consists of the normal vector fields on E. We define a metric on MT-L 2 {Q.) 
by means of 

<W(E 1; S 2 ) := d H {Af 2 Y>i,N 2 Y>2), 
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where dn denotes the Hausdorff metric on the compact subsets of R™ introduced 
in Section 2. This way MH 2 {Cl) becomes a Banach manifold of class C 2 . 

Let ds(x) denote the signed distance for £ as in Section 2. We may then define 
the level function ipy, by means of 

Ve(x) =0(d E (a!)), xeR n , 

where 

<H S ) = (!- X( s / a )) s gns + sx(s/a), sel. 
Then it is easy to see that £ = ^^(O), and V(^s(a;) = ^s(x), for x € S. Moreover, 
is an eigenvalue of V 2 (/3s(a;), and the remaining eigenvalues of V 2 < / 9s(a:) are the 
principal curvatures of £ at x £ £. 

If we consider the subset VWH 2 (r2,r) of MT-l 2 {Cl) which consists of all closed 
hypersurfaces F e M.% (0) such that T C satisfies a (interior and exterior) ball 
condition with fixed radius r > 0, then the map 

T : MH 2 (n 7 r) -> C* 2 (ft), T(T) := pr, (3.6) 

is an isomorphism of the metric space Ai% 2 {Vt,r) onto Y(A4'H 2 (r2, r)) C C 2 (f2). 

Let s — (n — l)/p > 2. Then we define 

W;(n,r) := {L G AW 2 (n,r) : Vr € IF p s (fi)}. (3.7) 

In this case the local charts for L can be chosen of class Wp as well. A subset 
A C Wp(n,r) is said to be (relatively) compact, if T(A) C Wp(fl) is (relatively) 
compact. 

As an ambient space for the state manifold SAA-y of the Stefan problem with 
surface tension we consider the product space C{G) x A4T-L 2 , due to continuity of 
temperature and curvature. 

We define the state manifolds SM 11 7 > 0, for the Stefan problem as 
follows. For 7 = we set 

SM := {(u,T) € <7(fi) x MW 2 : u £ W 2 ~ 2 l p {Sl \ F), F e IF p 4 " 3/p , (3.8) 
u > in fi, |#(u)] + cr^ = 0, J(u) ^ on T, [dft,u] e W£- a/p (r)}, 
and for 7 > 

SM^ := {(u, F) e C(Q) x MM 2 : u € IF p 2 - 2/p (0 \T),Te W^ 3/p , (3.9) 
u > in CI, (i(«) - [^(«)] - o^)([^>(«)] + o-W) = 7(u)[dft,«] on T}. 

Charts for these manifolds are obtained by the charts induced by M.% 2 {Q.), fol- 
lowed by a Hanzawa transformation. 

Applying Theorem 13. II or Theorem 13.21 respectively, and re-parameterizing the 
interface repeatedly, we see that yields a local semiflow on SM. 1 . 

Theorem 3.9. Let p > n + 2, er > and 7 > 0. Then problem 11. 1)) generates 
a local semiflow on the state manifold SM~f. Each solution (u,T) exists on a 
maximal time interval [0, t*), where f* = i*(uo,Lo). 
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Time weights. For later use we need an extension of the local existence results 
to spaces with time weights. For this purpose, given a UMD-Banach space Y and 
[i E (l/f>, 1], we define for J = (0,£o) 

KpJJ; Y) := {u G L Pt i oc (J; Y) : t^^u G K s p {J; Y)}, 

where s > and K G {H,W}. It has been shown in [31 that the operator d/dt 
in Lp^(J] Y) with domain 

D(d/dt) = oHljJ; Y)={ue H'JJ; Y) : u(0) = 0} 

is sectorial and admits an _ff°°-calculus with angle it/ 2. However, it does not 
generate a Co-semigroup, unless fi = 1. This is the main tool for extending the 
results for the linear problem, i.e. Theorems 13.31 and 13.51 to the time weighted 
setting, where the solution space E( J) is replaced by 

E M (J)=E |i) i(J) xE p , a (J), 

with 

E M ,i(J)= {v G H^(J;L p (ft)) n L p ^J;H 2 p {tl\V) : [v] = 0, = 0}, 

E^ j2 (J) : = ^p,{f _1 ^ 2p (>^; ip(^)) n %y 2p (^ Hp(T,))nL p ^(J; Wp _1 ^ p (E)),7 = 0, 
E M>2 (J) := W^ 1/2p (J; L P (E)) n L P>M (J; W^ 1/p (E)), 7 > 0. 
In a similar way, the space of data is defined by 
E M ,i(J) :=L p ^(J;L p (fl)), 

E^, 2 (J) := <^ 1/2p (J; L P (E)) n L PiP (J; W^-Vf (E)), 
F ft3 (J) := ^ 2 - 1/2p (J; £ P (E)) n L Ptfl (J; W^ 1 '^)), 
F„(J) := F M ,i(J) x F M , 2 (J) x F M , 3 (J). 
The trace spaces for u and p for p > 3 are then given by 

v e W^- 2 /p{n\Z), A e w; +2 "- 3 ^(S), Pl e M/ p 4 ^ 2 - 6 /p(E), (3.10) 

where for the last trace - which is of relevance only in case 7 = 0- we need in 
addition [i > 1/2 + 3/2p. Note that the embeddings 

E M ,i(J) C{J x ft) n C{J; C 1 ^)), E M>2 (J) C(J; C 3 (E)) 

require p > 1/2+ (n + 2)/2p, which is feasible since p > n + 2 by assumption. This 
restriction is needed for the estimation of the nonlinearities, i.e. Proposition 13.71 
remains valid for /i G (1/2 + (n + 2)/2p, 1). 

The assertions for the linear problem remain valid for such fj,, replacing E(J) 
by E M (J), F(J) by F M (J), for initial data subject to (13. 10[) . This relies on the 
fact mentioned above that d/dt admits a bounded _ff°°-calculus with angle 7r/2 
in the spaces L p ^(J;Y). Therefore the main results in Denk, Priiss and Zacher 
[5] remain valid for /j, G (1/p, 1). This has recently been established in [2"S] . As a 
consequence of these considerations we have the following result. 
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Corollary 3.10. Let p > n + 2, \x £ (1/2 + (n + 2)/2p, 1]. a > 0, and suppose 
that ip, 7 £ C 3 (0, oo), d £ C 2 (0, oo) such that 7 = or j(u) > 0, u £ (0, 00), and 

k(u) = — uip"(u) > 0, d(u) > 0, it € (0,oo). 

Assume the regularity conditions 

u e ^- 2 / p (r!\ro)nc(n), u >o, r e w* +2 ^ 3 / p , 

and the compatibility conditions 

(a) [^(uo)]+cW(r ) = 0, [d(u )cU ] e Wp 4 ^- 2 - 6/p (ro), as weH as the well- 
posedness condition Z(uo) 7^ ore To, m case 7 = 0. 

(b) (li>(u )l + an(T ))(l(u )-l^(uo)]-o-U(r )) =7(«o)[d(uo)^ ] incase 
7>0, 

TTien t/ie transformed problem (|2.1[) admits a unique solution z — (v,p) £ E^(0,r) 
/or some nontrivial time interval J — [0, r] . 77ie solution depends continuously on 
the data. For each 5 > i/ie solution belongs to E(5, r), i.e. regularizes instantly. 

4. Equilibria 

Suppose (u*,r») is an equilibrium for (|1.1[) . Then <9tit* = as well as V* = 0, 
and we obtain 

div(d(u*)Vu*) = in fl \ T* 

<9„ n it* = on dfl 

[it*] = on T* (4.1) 

l*l>(u*)] +aU{T if ) = on T* 

[d(u*)9j/U*] = on FV 

This yields u* = const, hence H(T*) — — [■0(u»)]/cr is constant as well. If T* is 
connected (and f2 is bounded) this implies that T* is a sphere Sr„ (xq) with radius 
i?* = er/[/0 («*)!■ Thus there is an (n + l)-parameter family of equilibria 

£ := {(u*,Sr,(x )) : u, > 0, < R* = ff/ty («•)], BrM C fi}. 

Otherwise, is the union of finitely many, say m, nonintersecting spheres of equal 
radius. It will be shown in the proof of Theorem 14 . 5 f vii) that £ is a C^-manifold 
of dimension (mn + 1) in W*(Q \ F*) x Wp~ 1/p {T*). 

4.1. Conservation of Energy. As we have just seen, the equilibria of Ijl.ljl 
are constant temperature, and the dispersed phase consists of finitely many non- 
intersecting balls with the same radius. To determine u and R, taking into account 
conservation of energy, we have to solve the system 

E(u,R) := IfixleiH + |fi 2 |e a («) + |E| - E , 

n — 1 

[V>(u)l + = 0. 

In order not to overburden the notation, we use (u, R) instead of («* , i?* ) . The 
constant Eq means the initial total energy in the system. Since % = —a/R we 
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may eliminate R by the second equation R = a/[ip(u)J, and we are left with a 
single equation for the temperature u: 

tp(u) := E(u,R(u)) = \Q\e 2 (u) - ^i?»[e( W )] + ^^R^\ u ) = ^, (4.2) 

n n — 1 

with ipo = Eq. Note that only the temperature range [^(u)] > is relevant due 
to the requirement R > 0, and with 

R m — sup{i? > : £1 contains m disjoint balls of radius R} 

we must also have R < R m , i.e. with h(u) — \ip(u)\ 

K u ) > 

With e(u) — ip(u) — uip'{u), i.e. = h(u) — uh'(u) 7 we may rewrite (f(u) as 

„(«) = |fl|e a (u) + c n + („ - l)n^) , 

where we have set c„ — m -. cr". 

n n(n — 1) 

Next, with 

ah'ju) h'(u)R\u) 

R M = ~ , 2 / \ = 

we obtain 

V '(u) = \n\e' 2 (u) - [e'C^HOxl +mu,„(-^ - ftl*)])^- 1 

= |fi|/t 2 (V) - [/C(«)]|fil| +771W„u/l'(u)i?"- 1 (u) J R'(u) 

.h'(u)R 2 (u) 



(K(u)\l) L2(n) -uh'(u)\n 



G 



au( K (u)\l) L2(n) ] l 2 (u)Ri(u)m 

P(u)R 2 ( U )\Z\ J (711 

with l(u) = uh'(u). It will turn out that in the case of connected phases the 
term in the parentheses determines whether an equilibrium is stable: it is stable 
if <p'(u) < and unstable if <p'(u) > 0; see Theorem 15.21 below. 

In general it is not a simple task to analyze the equation for the temperature 

h'(u) 



tp(u) 



unless more properties of the functions 62(11) and in particular of h(u) are known. 
A natural assumption is that h has exactly one positive zero u m > 0, the melting 
temperature. Therefore we look at two examples. 

Example 4.1. Suppose that the heat capacities are identical, i.e. [k] = 0. This 
implies 

«fc"(u)=«ty"(«)] =-[«(«)] = 0, 
which means that h(u) = h + h\U is linear. The melting temperature then is 
< u m — —ho/hi, hence we have two cases. 
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Case 1. ho < 0, hi > 0; this means l(u m ) > 0. 

Then the relevant temperature range is u > u m , as h is positive there. We assume 
now that e 2 is increasing and convex. As u — > m to + we have h(u) — > 0, hence 
93 (u) — > 00, and also £>(«) — > 00 for u — > 00 since e 2 (w) is increasing and convex. 
Further, we have 



93' (u) = |S7 je^ (w) — n(n — l)c„ 



(/to + ftiti)" 



+1 : 



which shows that <p(u) is strictly convex for u > u m . Thus ip(u) has a unique 
minimum uq > u m , <p(u) is decreasing for u m < u < uq and increasing for u > Uo- 
Thus there are precisely two equilibrium temperatures it+ G (110,00) and u~ G 
(u m ,uo) provided t^o > ¥?C"o) an d none if i^o < ¥>(«o)- The smaller temperature 
leads to stable equilibria while the larger to unstable ones. 

Case 2. /to > 0, hi < 0; this means l(u m ) < 0. 

Then the relevant temperature range is u < u m as h is positive there. As u — > u m — 
we have /t(u) — > 0+ hence tp(u) — > —00, and as u — > 0+ we have </?(«) — > <^(0) = 
|O|e 2 (0) + c„//tg _1 , assuming that £2(0) := lim M ^o+ £2(11) exists. Further, for u 
sufficiently close to zero <p'(u) is positive, since «2 = £2 > 0, and tp'(u) — > —00 
as u — > u m — . Therefore f'{u) admits at least one zero in (0,u m ). But there 
may be more than one unless £2(1*) is concave, so let us assume this. Let uq € 
(0,ti m ) denote the absolute maximum of <p{v) in (0,u m ). Then there is exactly 
one equilibrium temperature w* G (u , u m ) if <^o < and it is stable; there are 
exactly two equilibria u~ G (0,uo) and u+ G (uo,u m ) if <p(0) < ipo < ip(uo), the 
first one is unstable, the second is stable. If ipo > ip(uo), there are no equilibria. 

Note that in both cases these equilibrium temperatures give rise to equilibria only 
if the corresponding radius R(u) is smaller than R m . 

Example 4.2. Suppose that the internal energies q(u) are linearly increasing, i.e. 

€i{u) = di + KiU, i = 1,2, 

where m > 0, and now [k] ^ 0. The identity ej = ipi — uipl then leads to 

ipi(u) — cii + bill — Kilt In it, i = 1,2, 

where the constants i>j are arbitrary. This yields with a = [a] , /3 = [6] and (5 = [«] 

h(u) — a + flu — 5u\nu. 

Scaling the temperature by u = u w with fl — Slnu = and scaling h we may 
assume fl = and <5 = ±1. Then we have to investigate the equation f(w) = ipi, 
where 

<p(w) =cw + { ^j + (n - M») = ±(« + w Inw), 
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with c > and a, (pi € R. The requirement of existence of a melting temperature 
w m > 0, i.e. a zero of h(w) leads to the restriction a < 1/e. 

Actually, the requirement that the melting temperature is unique, i.e. that ft, 
has exactly one positive zero implies a < 0. Indeed, for a e (0, 1/e) there is a 
second zero w_ > of ft, and h is positive in (0,w_). Equilibrium temperatures 
in this range would not make sense physically. 

Also here we have to distinguish between two cases, namely that of a plus-sign 
where the relevant temperature range is w > w m and in case of a minus-sign it is 
(0,w m ). Note that h is convex in the first, and concave in the second case. 

Case 1: For the derivatives we get in the first case 



We have (p(w) — > oo for w — > oo and for to — > io m +, hence y(to) has a global 
minimum u in (u m ,oo). Further, ip"(w) > in (w m ,oo), hence the minimum is 
unique and there are precisely two equilibrium temperatures w~ £ (w m ,wo) and 
w+ e (w ,oo), provided ip\ > (p(w ), the first one is stable, the second unstable. 
To prove convexity of ip we write 



f(w) = 2w(ti(w)) 2 - h(w)(3 + ti(w)) = 2w(l +lnw) 2 - (a + wlnw)(4 + lnio). 
We then have f(w m ) = 2w m (l + lnw m ) 2 > 0, and 



for a < 1/e < w m < w. Let us illustrate the sign in h with the water-ice system, 
ignoring the density jump of water at freezing temperature. So suppose that 
consists of ice and fii of water. In this case we have K\ > K2 hence 6 < which 
implies the plus-sign for h. Here we obtain > w m , i.e. the ice is overheated. 
Equilibria only exist if 4>i is large enough, which means that there is enough energy 
in the system. If the energy in the system is very large then the stable equilibrium 
temperature w~ comes close to the melting temperature w m and then R(w) will 
become large, eventually larger than R*. This excludes equilibria in f2, the physical 
interpretation is that everything will eventually melt. On the other hand, if Oi 
consists of ice and of water, we have the minus sign, which we want to consider 
next. Here we expect under-cooling of the water-phase, existence of equilibria only 
for low values of energy, and if the energy in the system is too small everything 
will freeze. 

Case 2: Assume the minus-sign for h and let a < 0. Then the relevant 
temperature range is (0, w m ). Here we have tp(w) — > — oo as w — > w m — and 
(p(w) -> 1/H"- 1 > as w -)• 0+. 




(n + l)w(ti(w)) 2 - h(w)(3 + h'(w)) = (n - l)w{h'{w)f + f(w), 



where 



f'(w) = (1 + lnw;) 2 + 1 - a/w > 1 - a/w > 0, 



QUALITATIVE BEHAVIOR OF STEFAN PROBLEMS WITH SURFACE TENSION 25 



To investigate concavity of ip in the interval (0, w m ), we recompute the deriva- 
tives of ip. 

1 w(h'{w)) 2 - 



<p'(w) =c-(n-l){ — — + n / - /; \ 
y ' ' J lh n (w) h n+1 (w) J 



(w) h n+1 (w) 
V "{ W ) = n(n l)J^L{(n + l)w(h'(w)) 2 + h(w)(3 - h'(w))}. 

Setting w+ = 1/e, for w G (w + ,w m ) we have ft.(w) > and h'(w) < hence 
(p"{w) < 0. On the other hand, for w € (0, w+), both /i(w) and h'(w) are positive. 
Then we rewrite 

(n + l)w(/i'(w)) 2 + 3/i(w) - h(w)h'(w) = (n- l)w(l + lnw) 2 + f{w), 

where 

f(w) = 2w(h'(w)) 2 + h(w)(3 - ti(w)) 

= 2w(l + lnw) 2 — (a + w In w)(4 + lnw), 

f'(w) = (1 + lnw) 2 + 1 - a/w > 0, 

provided a < 0. This shows that / is increasing, f(w) — > — oo as m -> 0+, 
and /(1/e 3 ) = 11/e 3 — a > 0. On the other hand, the function w(l -I- lnw) 2 is 
increasing in (0, 1/e 3 ), hence p"(w) has a unique zero w_ € (0, 1/e 3 ). Therefore y> 
is concave in (0, w_)U (w + , w m ) and convex in (w_, w+), and <// has a minimum at 
w_ and a maximum at w + . Observe that tp'(w) < c, <p'(w) -oo for w ^ w m — 
and y>'(0) = c — (n — l)/|a|™ < tp'(w+). Therefore, <// may have no, one, two, 
or three zeros in (0,w m ) depending on the value of c > 0. However, if c > is 
large enough then ip' has only one zero w\ which lies in (w + , w m ). In this case <p 
is increasing in (0,wi) and decreasing in (wi,w m ), hence for ipi e (tp(0) , ip(wi)) 
there are precisely two equilibrium temperatures, the smaller leads to unstable, 
the larger to a stable equilibrium. If ip\ < ip{0) there is a unique equilibrium which 
is stable, and in case ip\ > <p(w\) there is none. However, in general there may be 
up to four equilibrium temperatures. 

4.2. Linearization at equilibria. The linearization at an equilibrium (u*,r*) 
with i?* = a/lip(u*)], reads 

K^dtv — d*Av = f in ft \ I\ 
d Va v = on <9f2 
{vj =0 on T* 
(h/u*)v + oA*p - 7*9 t p = g on T* 
l*d t p - [d*d u v} = h in 
v(0) = v , p(0) = p . 
Here 

1 n — 1 

k* = d* = d(u*), l*=l(u*), 7* = 7(u*), A* = -( — ^ — hA,), 



(4.3) 
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where A* denotes the Laplace- Beltrami operator on T*. 

We note that if Z» = and 7* = then the problem is not well-posed. On the 
other hand, if l t ^ and 7* = 0, then the operator — Lq defined by 

d(l ) = {(v, P ) g p p 2 (i!\r,)nc(ii)] x w^V'p*) : 

d va v = o, (i./u.)« + oA* P = o, [d.d v v] g w£- 2 /p(r*)}, (4.4) 

L (u,p) = ((-d*/n*)Av, — [(eJ*/i*)0„v]), 
generates an analytic Co-semigroup with maximal regularity in 

X Q := L„(fi) x H^ 2 - 2 / p (r»). 
More precisely, we have the following result. 

Theorem 4.3. Let 3 < p < 00, a > 0, suppose 7, = and Ze£ Z» =^ 0. Then for 
each finite interval J = [0,to], there is a unique solution z = (v,p) G E(J) of 
if and only if the data (/, g, h) and zq = (vq, po) satisfy 

(f, g ,h)e¥(j), z„e[if;-*(!i\r,)nc(!i)]xfl';-*(r,) 

and the compatibility conditions 

{h/u*)v + o-A*p = g(Q), h(0) + [d*d v v ] G M^ 2 - 6 /P(r t ). 

The operator —Lq defined above generates an analytic Co-semigroup in Xo with 
maximal regularity of type L p . 

In case 7, > 0, similar assertions are valid for L 1 in 

X 1 := L P (Q) x W^iT*), 

where 

d(l 7 ) = {( V ,p) e p p 2 (J!\r,)nc(ii)] x w£-Vp(r.) : 

9 1 ^w = 0, (/ 2 /it*)« + l*aA*p = 7*[d*^u]}, (4.5) 
L 7 (v,p) = ((-d*/K*)Av,-(a/'y*)A tt p - (l*/u*j*)v). 

The main result on the problem (|4.3|) for 7* > is the following. 

Theorem 4.4. Let 3 < p < oo, and suppose a, 7* > 0. TTien /or eac/i finite 
interval J — [0, to],. Ace is a unique solution z — (v,p) G E(J) 0/ i/ and 

onZy if t/ie data (/, g, h) and z Q = (vq, po) satisfy 

(/) 9i h) g f(j), z,e[if3\r,)nc(!i)]xH';-*(r t ) 

and t/ie compatibility condition 

(ll/u*)v + l*aA*p - 7*[d9 l/ w ] = l*g(0) + 7*/i(0). 

77ie operator —L 1 defined above generates an analytic Co-semigroup in Xj with 
maximal regularity of type L p . 
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Proof of Theorem \4-3\ and Theorem \4-4\ These results are, up to the last asser- 
tions, special cases of Theorems 13.31 resp. 13.51 in Section 3. In addition, since the 
Cauchy problem for L 7 has maximal L p -regularity, we conclude in both cases by 
[29l Proposition 1.2] that — L 1 generates an analytic Co-semigroup in X 1 . Recall 
that the spaces E( J) are different for 7 = and 7 > 0. □ 

4.3. The eigenvalue problem. By compact embedding, the spectrum of L 1 
consists only of countably many discrete eigenvalues of finite multiplicity and is 
independent of p. The eigenvalue problem reads as follows 

k*\v — d*Av = in f2 \ T* 

d Vu v = on dfl 

[v] = on (4.6) 

(l*/u*)v + oA^p — 7*Ap = on T* 

l*Ap — ld*d v vj —0 on r*. 

Assume first that is connected. As shown in [32], A = is always an eigenvalue, 
and N(L 7 ) is independent of 7* > 0, k* and d*. We have 

N(L 7 ) = span |(£^, -1), (0, Fi), . . . , (0, F„)| , (4.7) 

where the functions Yj denote the spherical harmonics of degree one, normalized by 
{Yj\Yk) l 2 (t ,) — $jk- N(L 7 ) is isomorphic to the tangent space of £ at G £ . 

Let A 7^ be an eigenvalue with eigenfunction (v,p) 7^ 0. Then (|4.6|) yields 

\{\^v\l 2(n) - au*(A*p\p) L3( rj} + \\fd~^J v\\ 2(n) + 7*"*|A| 2 |Hi 2 (r,») = 0. 

Since A* is selfadjoint in L2(r*), this identity shows that all eigenvalues of L 7 are 
real. Decomposing v = v Q + v, p = p + p, with (K*|wo)i 2 (n) = (po|l)i 2 (r,) = °j 
this identity can be rewritten as 



A{|v^w || 2 (n) - t7U*(A*pti\po)L a {r.) + A "*7*|PoL 2 (r»)} + I V d *^ v o\l 2 (n) 
+ [A7*w» + Z*|r*|/(/c*|l) i2 (n) — o-u*/Rl]\p 2 \rJ — 0. 

In the case that r„ is connected, the bracket determines whether there is a positive 
eigenvalue. 

If T* = Ui<i< m I* consists of m > 1 spheres r$, of equal radius, then 

Ar(i 7 )= sp an|(£^,-l) ) (0,y 1 i ),...,(0,^): 1 < I < m| , (4.8) 

where the functions Yj denote the spherical harmonics of degree one on rj, (and 
Yj e on LL^rj.)) normalized by {Yj\Y£) La (pij = 8 ]k . N(L y ) is isomorphic to 
the tangent space of £ at (u*,r*) e £, as will be shown in Theorem 14 . 5 1 below . 
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Theorem 4.5. Let a > 0, 7* > 0, h 7^ 0, and assume that the interface T* 
consists of m > 1 components. Let 



(ii) is a an eigenvalue of—L 7 with geometric multiplicity (ran + 1). 

(iii) is semi-simple if 7^ 1 . 

(iv) If r* is connected and < 1, then all eigenvalues of —L 7 are negative, 
except for 0. 

(v) // C* > I? then there are precisely m positive eigenvalues of —L 1 , 

(vi) If < 1 iften — L 7 ftas precisely m — 1 positive eigenvalues. 

(vii) N(Lj) is isomorphic to the tangent space Tt Ut ,T„)£ of £ at (u*,r«) G £. 

Remarks 4.6. (a) The result is also true if I* — and 7* 7^ 0. In this case 
ip'(u*) = (k + I l)z, 2 (n) > and = 00, hence the equilibrium is always unstable. 

(b) Note that does neither depend on d, nor on the undercooling coefficient 7. 

(c) For the Mullins-Sekerka problem, that is, for k = 0, we have = 0, in 
accordance with the result obtained in |13j . 

(d) It is shown in [32 that in case = 1 and connected, the eigenvalue is 
no longer semi-simple: its algebraic multiplicity rises by 1. This is also true if T* 
is disconnected. 

Proof of Theorem \4-5\ For the case that T* is connected this result is proved in 
[32] , The assertions (i)-(iii) also remain valid in the disconnected case. However, 
the proof of [32l Theorem 2.1(e)], addressing instability, is not completely correct, 
as it relies on the assertions [521 Proposition 3.2(b) and Proposition 5.1(c)] which 
are incorrect. (We remark, though, that the instability result of [3SJ Theorem 1.3] 
indeed is valid.) Here we give a modified proof for ;32] Theorem 2.1(e)] which also 
applies in case T* is not connected. 

It thus remains to prove the assertions in (v), (vi), and (vii). If the stability 
condition f» < 1 does not hold or if T* is disconnected, then there is always a 
positive eigenvalue. To prove this we proceed as follows. Suppose A > is an 
eigenvalue, and that p is known; solve the heat equation 



g-M«(AC«|l).L 2 (n) 

^|r»| 



(4.9) 




k*\v — d*Av = in f2 \ T* 



d Vn v = on dfl 
[v] =0 on r» 
— = h on 



(4.10) 



to get v = S\h, with S\ being the solution operator. Then taking the trace at 
r* we obtain v\r, = N\h, where N\ denotes the Neumann-to-Dirichlet operator 
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for the heat equation (|4.10p . Setting h — —\l*p this implies with the linearized 
Gibbs-Thomson law the equation 

[{l 2 Ju^XN x +^X]p-aA,p = Q. (4.11) 

A > is an eigenvalue of — L 1 if and only if (|4.11j) admits a nontrivial solution. 
We consider this problem in L2(T*). Then A* is selfadjoint and 

-<j{A*g\g) L ^ v ,) > -^2bl! 2 (r»)- 

On the other hand, we will see below that N x is selfadjoint and positive semi- 
definite on ^(r*). Moreover, since A* has compact resolvent, the operator B x := 
l(l%/u*)\N\ +7»A] — a A* has compact resolvent as well, for each A > 0. Therefore 
the spectrum of B x consists only of eigenvalues which, in addition, are real. We 
intend to prove that in case either T* is disconnected or the stability condition 
does not hold, B Xo has as an eigenvalue, for some Ao > 0. 

We will need the following result on the Neumann-to-Dirichlet operator N x . 
We denote by e the function which is identically to one on T*. 

Proposition 4.7. The Neumann-to-Dirichlet operator N\ for problem (14. 10)) has 

the following properties in ^(T*)- 

(i) If v denotes the solution of (|4.10p . then 

(N x h\h) L2(rt) = A|VML(n) + lv^V«|i s(n)) A > 0, h £ L 2 (T*). 

(ii) For each a £ (0, 1/2) and Ao > there is a constant C > such that 

(N x h\h) L2{rr) > ^\N x h\l 2(Tr) , h £ L 2 (T„), A > A . 
In particular, N x is injective, and 

C 

\ N \\i3(L 2 (r,)) < -j—, A > A . 

(iii) On L 2: o(r*) = {0 £ £a(r*) : (0\e)L 2 {r,) = 0}, we even have 

(N x h\h) L2{Tr) > (1 ^ A I" \N x h\l 2(Tt) , h £ L 2 ,o(r,), A > 0, 
and 

C 

\N\\B(L 2 , (r,),L 2 (T,)) < 7j - A ^ Q > A > 0. 

In particular, for A = 0, (|4.10p is solvable if and only if (/i|e)r„ = 0, and then the 
solution is unique up to a constant. 

Proof of Proposition \4- 7[ The first assertion follows from the divergence theorem. 
The second and third assertions are consequences of trace theory, combined with 
Poincare's inequality. The last assertion is a standard statement in the theory of 
elliptic transmission problems. We refer to [35]. □ 
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Proof of Theorem \4-5\ continued: 

(a) Suppose first that T* is connected. Consider h — e. Then with c* :=l1/u* > 
we have 

(S A e|e) i2(r ,) = c*A(iV A e|e) L2(IM + |e|| 2(r><) - |e|| 2(r>) . 
We compute the limit lim A _^ A(iV A e|e)x, 2 (r < ) as follows. First solve the problem 



d vn v = 



d*Av — — K*ao in O \ 

on 3S1 

on r* 

on r*, 



M = o 

. -{d*d v v] = e 



(4.12) 



where ao = |r*|/(K*|l)i 2 (n) ; which is solvable since the necessary compatibility 
condition holds. Let vq denote the solution which satisfies the normalization con- 
dition (k* |vo)l 2 (si) = 0. Then v\ := S\e — vq — ciq/X satisfies the problem 



' k*\v\ — d*Av\ = — k^Xvq in ft \ r* 

d vn v = on dfl 

[v x ] = on T* 

-[d*d u V\} =0 on r*. 



(4.13) 



By the normalization (k* |i>o)l 2 (o) = we see that v\ is bounded in W% (fi \ r*) 
as A — > 0. Hence we have 

lim XN\e = lim[Au A |r» + Au |r» + do] = ao = |r*|/(«;*|l)z, 2 (n)- 

A^O A-s-0 



This then implies 



ir. 



^|r*l<o, 



lim(B A e|e) L2(rt) = c, . ■ 

l K *|J-jL 2 (n) 

if the stability condition does not hold, i.e. if C* > 1< 

(b) Next suppose that is disconnected. If F* consists of m components T^, 
fc = 1, m, we set e& = 1 on and zero elsewhere. Let h = ^ fe a^ek ^ with 
X) fe ak = 0, hence Qo^i = h, where Qo is the canonical projection onto L2,o(r*), 



Qoh = h — 



(h\e) 



L 2 (r, 



Then 



lim \N\h = lim \N\Q Q h = 0, 

A^0 A^O 



since N\Q is bounded as A — > 0. This implies 
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(c) Next we consider the behavior of (B x h\h)i J2 ^r t ) as A — > oo. We want to show 
that B\ is positive semi-definite for large A. For this purpose we introduce the 
projections P and Q by 

m 

Ph = c m ^2(h\e k ) L2{rr) e k , Q = I - P, 

fe=i 

where c m = m/|r*| in case T* has m components. Recall that |r*| = |r*|/m for 
k = l,...,m. Then with h k = (h\e k )L 2 (r,) 

\(N x Ph\Qh) L2irt) \ < c m J2\h k \\(N>Qh\e k ) L2{rt) \ 

k 

< CY,\hk\\N x Qh\ L2(rm) < C\- a ' 2 \hk\(NxQh\Qhf^ {Tt) 

k k 

< CX^lY, \hk\ 2 + m(N x Qh\Qh) L2irt) } 

k 

< C\- a / 2 [\Ph\l 2{Vt) + (N x Qh\Qh) L2{r ,)}, 

for A > 0, and C standing for a generic positive constant, which may change from 
line to line. Hence for A > Ao, with Ao sufficiently large, we have 

(N x h\h) L2{rr) = (N x Qh\Qh) L2(rt) +2(N x Qh\Ph) L2{rr) + (N x Ph\Ph) ia(r .) 

> \(N x Qh\Qh) L2[Tr) + (N x Ph\Ph) L2(Tr) - JL\Ph\l 2(rt) . 

This implies 



(B x h\h) L2( p r) =c*\{N x h\h) L2(rt) +7*X\h\l 2(Tr) - (A„h\h) L ^ r ,) 
~2 



> ^{NxQh\Qh) L2(r ^ + c*\(N x Ph\Ph) L2iTt) 



-(A*Qh\Qh) L2( r m )-c\Ph\l 2(r ^. 

Since N x is positive semi-definite and also —A*Q has this property, we only need 
to prove X(N x Ph\Ph) L2 (^ rt - j — > oo as A — > oo. 
To prove this, similarly as before we estimate 

KiV^le,)^,)! < C|A^ A e,U 2(r ^ < CX^ 2 {N x e^ 2(TtV 
and choosing Ao sufficiently large this yields 

r C i 

(N x Pg\Pg) L2{rt) > c min(JV A e i je i ) i3(r . ) - — ^ |^slL(r,v 



,a/2 
A 



Therefore it is sufficient to show 



lim X(N x e k \e k ) L2 , r , = oo, k=l,...,m. (4.14) 

So suppose, on the contrary, that Xj(N x g\g)i J2 (r t ,) is bounded, for some g = e k and 
some sequence Xj —¥ oo. Then the corresponding solution Vj of (|4.10l) is such that 
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XjVj is bounded in Z>2(^), hence has a weakly convergent subsequence. W.l.o.g. 
XjVj — > Voo weakly in ^(Q). Fix a test function ip G T>(Q \ T*). Then 

^( k *^#)l 2 (o) = {d^Avjlip)^^ = (w J |d„A?/')L 2 (n) 

= (XjVjldtA^L^/Xj 

as j — > oo, hence = in L>2{Vt). On the other hand we have 



0<|r»|/m = J gdT* = J -ldd uVj ]dr 



dAvjdx = Xj / n^Vjdx — ¥ / n^Voodx, 
in Jn Jn 

hence is nontrivial, a contradiction. This implies that (|4.14[) is valid, provided 

h > 0. 

On the other hand, in case I* — we have 7, > 0, hence A-^* l/?lL 2 (r ) ~~ ^ 00 > so 
also in this case B\ is positive semi-definite for large A. 

(d) Summarizing, we have shown that B\ is not positive semi-definite for small 
A > if either r» is not connected or the stability condition does not hold, and 
B\ is always positive semi-definite for large A. Set 

Ao = sup{A > : is not positive semi-definite for each /1 £ (0, A]}. 

Since B\ has compact resolvent, B\ has a negative eigenvalue for each A < Ao. 
This implies that is an eigenvalue of B\ 01 thereby proving that — L 7 admits the 
positive eigenvalue Ao- 

Moreover, we have also shown that 

Boh = lim c*XN\h — aA*h — cAT A I L tQ)P h — aA*h, 

where PqH := (I — Qo)h = (/i|e)i 2 (r,)/|r*|. Therefore, B has the eigenvalue 
c*|r* |/(k*|1)l 2 (s7) — a /R* with eigenfunction e, and in case m > 1 it also possesses 
the eigenvalue — cr/i? 2 with precisely m — 1 linearly independent eigenfunctions of 
the form a^k with X)fc a fc = 0- This implies that — L 1 has exactly m positive 
eigenvalues if the stability condition does not hold, and m — 1 otherwise. 

(e) It remains to show assertion (vii). Suppose for the moment that T* consists 
of a single sphere of radius i?* = a/ltp(u*)}, centered at the origin of R". Suppose 
S C SI is a sphere that is sufficiently close to T*. Denote by (z\, . . . , z n ) the 
coordinates of its center and let zq be such that <r /[■;/>(«* + zo)] corresponds to its 
radius. Then, by [131 Section 6], the sphere S can be parameterized over T, by 
the distance function 



P( z ) = z i Y i ~ R * 

3=1 



^ i=i i=i 
Denoting by O a sufficiently small neighborhood of in the mapping 

[z ^ *(z) := (u, + z , p{z))\ : O -► % 2 (ft) x W^"^^) 
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is C 1 (in fact C k if tp is C fc ), and the derivative at is given by 

n 

*'(0)/l=(l,-<7[^K)]/^(M*)] 2 )/lO+(0,^^) ! ^ K " +1 

3=1 

Noting that £t[V>'(u*)J/I[V'( u *)] 2 = Z*i? 2 /(<7U*) we can conclude that near 
the set £ of equilibria is a C 1 -manifold in Wp(Cl) x Wp (r*) of dimension n+1, 
and that the tangent space T( u ,,r»)(^) coincides with N(Lj), see (|4.7[) . It is now 
easy to see that this result remains valid for the case of m spheres of the same 
radius i?*. The dimension of £ is then given by (inn +1), as mn parameters are 
needed to locate the respective centers, and one additional parameter is needed to 
track the temperature. □ 



5. Nonlinear Stability and Instability of Equilibria 



Before we discuss nonlinear stability of equilibria, we need some preparations. 
The first observation is that the equations near an equilibrium (u* , T» ) € £ can 
be restated as 



K*d t v — d^Av = F*(v, p) in fi \ 1% 

H =0 on r* 

(l*/u*)v + aA*p — J*dtp = G*(w, p) on T* 

Udtp— [d* = H t (v,p) on r» 
u(0) = w , jo(0) = p . 



(5.1) 



where 



F*(u, p) — (k» — k(u» + v))d t v + (d(u* + v) — d*)Av + d(u* + v)M 2 {p) : V 2 v 
- <f («, + v)\(I - Mi(p))Vv\ 2 + d(u, + u)(M 3 (p)|Vu) 
+ k(u* + v)ll(p)(u* + v), 
([V>(tt* + u)] + crH(p)) + {l*/u*)v + aA*p+ (j(u* + v)f3(p)-j*)d t p, 
m* + u) — d*)<9„v] + (/» — Z(m» + v))d t p 
u* + v)Vv]|M 4 (p)V s p) + 7 (u. + v)/3(p)(<9 t p) 2 , 



G*{v,p) 



see Sections 2 and 3 for the definition of Mj(p), j = 1, ... ,4. The nonlineari- 
ties satisfy F„(0,0) = G«(0,0) = if, (0,0) = as well as i^(Q,0) = G'„(0,0) = 
^(0,0) = 0. 
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The state manifold for problem (|5.1J near the equilibrium (it*, T*) can then be 
described by 

sm = {(v lP ) e w^p(n\r,) x w^'n?*) ■ d vn v = o, H = o, (5.2) 

(!,/u,)v+trA*p = G.(v,p), [<Udvv]+H*(v,p) G W^ 6/p (T,)}, 
for 7* = 0, in case /» ^ (otherwise the linear problem is not well-posed) and 

sm^ = {(v, P ) e w*- 2 / p (n \ r.) x % 4 " 3 / p (r,) : 9„ n <; = o, [«] = o, (5.3) 

(ll/u*)v + l*aA*p - 7*[rf*9iyt'l = l*G*(v,p) + j*H*(v,p)}, 

in case 7* > 0. 

We would like to parameterize these manifolds over their tangent spaces at (0, 0), 
given by 

z = {(v,p) G [w;-^(!]\r,)nc(si)] x W^(T t ) : (5.4) 

d Vn v = 0, + oA*p = 0, G W 2 ~ 6/p {T,)} , 

respectively, for 7* > 

z 1 = {(s,p) g [VK p 2 - 2 /p(fi\r,)nc(n)] x w p 4 - 3 /p(r«) : (5.5) 

d„ n i; = 0, {il/u*)v + UaA*p - 7*[d*d„t;] = 0}. 
Note that the norm in Z 1 for 7 = is given by 

\(v,p)\z = \v\ w 2-*/p + \p~\ w *-Vp + 
while for 7 > it is given by \{v,p)\^ = \v\ w 2-2/ P + \p\ w 4-?,/ P . 

In order to determine a parameterization, we consider the linear problem 

k*u!v — d*Av = in f2 \ r* 
= on <9£1 
H =0 on r* (5.6) 
(U/u*)v + <tA*p — j*wp = g on 

l*u)p—\d*d v v\=h on r*. 

We have the following result. 

Proposition 5.1. Suppose p > 3, 7* > 0, /* 7^ in case 7* = 0, and lj > 

is sufficiently large. Then problem ()5.6|) admits a unique solution (v, p) with 
regularity 

v g w^- 2 /p(n XT,), P e % 4 - 3 /p(r,) 

i/ and only if the data (g, h) satisfy 

The solution map [(g, h) 1— >• (u, p)] is continuous in the corresponding spaces. 
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Proof. This purely elliptic problem can be solved in the same way as the corre- 
sponding linear parabolic problems, cf. Theorems 4.3 and 4.4. □ 

For the parametrization we pick u) > sufficiently large. Given z — (v, p) £ Z^ 
sufficiently small, we can solve the auxiliary problem 



k*ujv — d,Au 






M = o 

(i*/it*)v + aA*p — j*ujp = G*(v + v, p + p) 
Uuip — \d»d v v\ = H*(v + v, p + p) 



m 
on 
on 
on 
on 



on 

r, 
r* 
r. 



(5.7) 



by means of the implicit function theorem, employing Proposition (5TTJ This yields 
a unique solution z = (y, p) = <j>(z) G Wp~ 2 (Cl \ r») x Wp~~ 3 ^ p (T*) with a C 1 - 
function 4> such that 0(0) = as well as ^'(0) = 0. One readily verifies that 
z = z + 4>(z) £ SM-y. To prove surjectivity of this map, given (v,p) £ SM-y, we 
solve the linear problem 



K*U1V - 



d*Av 






M = o 

(h/u*)v + oA*p — -f*ujp — G*(v, p) 
l*up- \d*d v v\ = H*(v,p) 



in 

on 
on 
on 
on 



n\r, 
on 
r* 
i 

r* 



(5.8) 



and set z 



Then z £ Z~ and z 



hence the map [z4z + 4>{z)} is 



also surjective near 0. We have, thus, obtained a local parametrization of SM. 1 
near zero over the tangent space Z 1 . 



Next we derive a similar decomposition for the solutions of problem (|5.ip . Let 
z o = (So)<K^o)) e SM 1 be given and let z £ E(a) be the solution of (15.1[) with 
initial value z . Then we would like to devise a decomposition z — z oc + z + z, 
where z(t) £ Z 7 for all t £ [0, a], and where = Zoo + 4>(Zoo) is an equilibrium 
for (|5.1I) . In order to achieve this, we consider the coupled systems of equations 

k*ujv + K*d t v — d*Av = F tf (z ao + z + z) — F^(z OQ ) 



d Un v = 
[v]=0 

(l*/u*)v + o-A*p - 7*(d t p + ujp) = G*(zoo + z + z) - G»(z oc ) 
hup + hd t p - {d*d v v\ = H*(zoo + z + z) - H„(z 00 ) 

z(0) = <p(z Q ) - ^(Zoo), 



(5.9) 



and 
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(5.10) 



It should be mentioned that P*(Zoo) = 0, as can be seen from the equilibrium 
equation for (|5.ip and the fact that = constant for = (woo,Poo)- For 
reasons of symmetry and consistency we will, nevertheless, include this term. 
Equations (|5.9[) (|5.10|) can be rewritten in the more condensed form 



where we use the abbreviation L 7iW to denote the linear operator on the left hand 
side of (|5.9[) , and N to denote the nonlinearities on the right hand side of (|5.9I) . 
respectively. We are now ready to formulate the main theorem of this section. 

Theorem 5.2. Suppose a > 0, 7* = 7(1/*) > and h — l(u*) ^ in case 7* = 0. 
Then in the topology of the state manifold SAi-y we have: 

(a) (tt*,r*) £ £ is stable i/T* is connected and < 1. 

Any solution starting in a neighborhood of such a stable equilibrium exists 
globally and converges to another stable equilibrium exponentially fast. 

(b) (u*,r*) £ £ is unstable ifT* is disconnected or if > 1. 

Any solution starting and staying in a neighborhood of such an unstable 
equilibrium converges to another unstable equilibrium exponentially fast. 

Proof, (a) We begin with the case that (it*, T*) is linearly stable. Then according 
to Theorem l4.5l we have = N(L 7 ) © R(L 7 ). Let P c denote the projection onto 
X^ :— iV(L 7 ) along X^ :— R(L~,) and P s = I — P c the complementary projection 
onto J?(L 7 ). We parameterize the set of equilibria £ near over JV(L 7 ) via the C 1 - 
map [xh+x + ip(x) + <j>(x+ip(x))] such that ip(0) = ip'(0) = and 0(0) = <f>'(0) = 0. 
It follows from the equilibrium equation associated to (|5.1[) (recall that F*(z e ) 
vanishes for any equilibrium z e ), and from the definition of <f> that the mapping ip 
is determined by the equation 



Since L 7 is invertible on X*, ip 6 C 1 (Bx^(r), D(L^)) is well-defined by the implicit 
function theorem and ip(0) = ip'(0) = 0. 

For Xqo € X^ sufficiently small we set z^ := Xoo + ip( x oc) + 4>(xoc + ?/>(x<x>))- 
Then z^ is an equilibrium for (|5.1[) and we will now consider the decomposition 



h 7tU1 Z = N(Zoo+ Z + z) - N(Zoo), 

Z + LjZ = UJZ, 



Z(0) = <j)(z ) ~ 4>{5oo), 
Z(0) = Z - Zoo, 



(5.11) 



L^(x)=P s u<f>(x + ip(x)), xeB x .(r). 



(5.12) 
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z = Zoo + z + z introduced in (|5.9|) - (|5.10|) . or (|5.11|) . respectively. With the ansatz 

z = x + ^(x 00 + x)-V(x o) + y, (x,y)ex;xi; (5.13) 

for x, Xqo g X^ small enough, the second line in (I5.11|) becomes 
j x = P c ujz, x(0)=x -x oo , 

|y + i;y = 5(x 00 ,x,z), y(0) = y , l "" Ul 



where 



and 



S^XoojX, Z) = P S (JZ - ^'(Xoo + x)P c uz - L*[V>(Xoo +x) - ^(Xoo)], 



z = x + ^(x ) + y , (xo.yjlcI^x^nZj. (5.15) 

Next we show that the system of equations (|5.14|l admits a unique global solution 
(xqo , x, y) , where the functions (x, y) are exponentially decaying, provided z is ex- 
ponentially decaying and (xo, yo) is sufficiently small. For this let us first introduce 
some more notation. For S > we set 

Ei(R + , S) := {v : e 5t v € E, (R+)}, i = 1,2, 

F i (R +J <y) := {« : e 5t v G F,(R+)}, j = 1,2,3, 

endowed with the norms 

Mk(R + ,S) = l|e <5 *^|| Ei (M + ), |w|f 3 (r + ,5) = |e 5t f|F J (R + )- 

The spaces E(R+,d) and F(R + ,<5) are then defined analogously as in Section 3. 
We also need the space 

X(R+,c5) := iJ p 1 (R + ,<5; X 7 ) n L P (R + ,S; D(L y )), (5.16) 

where H£(R + ,6; E) denotes all functions v : R + — > E such that e Si v € H£(M. + ;E), 
with E a given Banach space. Finally, let 

M 1 (r,5): = {(x ,y ,z)eXc x {X* n Z 7 ) x E(R+, <5) : |(x ,y )|^ < r}. 

For given (xq, yo, z) G Bi(ro, S), with rp sufficiently small, we set 



x + / P c ujz(t) dr, 
o 

P c ujz(r)dT, (5.17) 

y := f^ + L 7' tr ) (<5(xoc,x,z),y ). 

Here we used the notation trio := w(0). It should be observed that the functions 

(xoo , x) occurring in the third line of (|5.17|) are defined through the first two lines 
in (|5.17j) . We now set 

&(x ,yo,z) := (xoo,x,y), (xo, yo, z) G Bi(Vo, S), (5.18) 

where tq is chosen sufficiently small. 
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Next we will show that there exists a number So > such that for any 5 £ [0, So] 
the mapping 6 has the following properties: 

6 G C(Bi(r ,<5),X^ x X c (R+,<5) x X s (R+,<5)), 6(0) = 0, (5.19) 

where 

X c (R+,<5) :=Hl(R + ,S;X°), X s (R+,<5) := <5; X») n L p (R+,S; D(L S 7 )). 

Writing 6 = (©1,62)63) we readily observe that 

61 eC°°(li(ro,i),^), 6i(0) =0. (5.20) 
For g £ L p (R + ,S;Xc), let (Kg)(t) := J f °° g(r) dr and note that 

/OO 
e ^)e^(r)dr. 

Young's inequality for convolution integrals readily yields 

K £ B(L P (R +7 S; X°), H^(R + ,S; X°)) , 

and this shows that 62 £ X C (R + , S). Hence we have 

6 2 G C oo (M 1 (r ,S),X c (R + ,S)) ) 6 2 (0) = 0. (5.21) 

Concerning the function 63, we know from Theorem l4.5f v) that s(—Lt), the spec- 
tral bound of (— Lt), is negative. Fixing Sq > with s(—Lt) < — So it follows from 



semigroup theory and the L p -maximal regularity results stated in Theorem 14.31 
and Theorem 14.41 that 



(j t +L s 1 ,tr)- 1 £B{L p (R +7 S;X^)xZ^X s (R +7 S)), Se[0,S o ], (5.22) 

where Zt = X* n Z 7 . This in conjunction with (|5.20|) - (|5.21|l and the definition of 
S implies 

6 3 £C(M 1 (r Q ,S),X s (R + ,5)), 6 3 (0) = 0. (5.23) 
Combining (f5T20]) - (j5T23j) then yields (j5~T9j) . 

For given (xo,yo,2) £ Mi(r ,S) let (x m ,x, y) = 6(x ,yo,z). Then we have 



x(t) 



- / P c ujz(t) dr = - P c ujz(t) dr + P c ujz(t) dr 
Jt Jo Jo 

x - + / P c ujz(t) dr, 
Jo 



thus showing that x solves the first equation in f|5 . 14[) . In summary, we have shown 
that (Xoo,x, y) = 6(x ,yo,z) is for every (x ,yo,z) £ Bi(r ,<5) the unique solution 
of in X* x X c (R + ,<5) x X s (R + ,<5), where S £ [0,S Q ]. 

Setting 

z = 3(x , yo, z) := x + ?A(xoo + x) - ^(xoo) + y, 

Zoo = 3oo(x ,yo,^) := Xoo + ?/K X oo) + 0(Xoo + ^(Xoo)) 
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for (xoo,x, y) = ©(xo,yo,z), we see that 

3eC(B 1 (r ,<T) ) X(M + ,<5)), 3(0) = 0, 

and 

3oa eC(Bi(ra, 5), Zoo), 3oo(0) = 0, (5.25) 

where Z x = [W|(ft\r»)n(7(ft)] x Wp~ 1/p (T*). It then follows from the derivation 
of f5TT3]l - (^Tgjl that 

(zoo,z) = (3(x ,yo,z),3(xo,yo,z)) 

is for every given (xo,yo,z) G V>i(r ,5) the unique (global) solution of (|5.10l) with 
z in the regularity class X(R + ,<5). In a next step we shall show that z in fact has 
better regularity properties, namely 

3eC(Bi(r 0) (S),E(R + ,<J)) ) 3(0) = 0. (5.26) 

In order to see this, let us first consider the case 7 = (which implies 7, = 0). 
From the fourth line of (I5.10p . the fact that z e X(R + , <5), and 

[v 1 ^ «| r .] G B(X(R + J),W^ 1/2p (R+,5; L P (T*))) 

follows 

p = (fi-aA^dh/u^v + fxp) e VK p 1 - 1 / 2f, (R + ,<5;i/ p 2 (r»)), 

where p, is in the resolvent set of a A* . From the fifth line of (|5 . 10|) , the fact that 
(z, z) G X(R + , 5) x E(R +) S), and trace theory for v follows 

implying that p e Wp /2 ~ 1/2p (R+, 5; L P (T*)). Hence holds for 7 = 0. 

If 7 > (and thus 7, > 0) we use the embedding 

iJ p 1 (R + , ( 5;W^- 1 /P(r,))nL p (R + , ( 5;^ p 4 - 1 /P(rO) I^- 1/2p (R + , <5; ff 2 (r*)) 

and the fourth equation in (|5.10[) to conclude that p € Wp 2 _1 ^ 2p (R + , (5; i p (r»)). 
Hence (|5.26|) holds in this case as well. 

Let us now turn our attention to equation (|5.9p , or equivalently, the first line of 
(|5.1ip . In a similar way as in the proof of }25] Proposition 10] (extra consideration 
is needed in order to deal with the additional terms involving Zoo ) one verifies that 
the mapping 

[(*«,, z) 1 y N( Zoo + z) ~ N( Zoo )] : U(<5) -> F(R+, 5) 

is C 1 and vanishes together with its Frechet derivative at (0, 0). Here V(S) denotes 
an open neighborhood of (0, 0) in x E(R+, S). Let 

M(r,S) ={(xq,yq,z) e x {X" n Z-f) x E(R+,<5) : |(xo,yo,^)|[2,]2xE(K + ,i5) < ''0}, 

and let ext 5 <= B(Wp~ 2/p (Q \ r») n C(Q)) x Wp 4 ~ 3/p (r*), E(R+, <5)) be an appro- 
priate extension operator with (extan;o)(0) = wq. 
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For (xo,yo,z) G M(ro,S), with ro sufficiently small, we define 

Af(xo, yo, 2) := ^(^oo + 2 + ext 5 (^(2 ) - 0(*») - 5(0)) + z) - AT(*«>)- 

It follows from (IOI)) - (fOa that M G C(B(r , <5), F(R+, <5)), M(Q,0,0) = 0, and 
D 3 M(0,0,0) = 0. Moreover, 

M(x ,y ,z)(0) = N(zo)-N(z ao ), (x , y , z) e B(r , <5), (5.27) 

where we recall that z(0) = zq — Zoo, zo = zo + </>(zo), and Zoo = Zoo + <f>(zoo)- 
Finally, for (xo,yo,z) G B(ro,<5) let 

K(xo,y ,z) := (L 7iW , tr) -1 (M(xo, y , z), 0(z o ) - (5.28) 

It follows from (|5.27|l and the definition of <f> that the functions 

(M(x , y , z), 4>{z Q ) - (pizoo)) 

satisfy the necessary compatibility conditions, whenever (xo,yo,z) G B(ro,J). 
Slight modifications of the results in [9] then imply that K : B(ro,<5) — > E(R+,<$) 
is well-defined, provided cj is large enough (and S is in [0, Jo] w hh So as above). 
From the properties of the mappings N, ip and <fi, the definition of So and Zoo 
(recall that zo = xo + ?/>(xo) +yo, Zoo = Xoo + V'( X °°))> and the contraction mapping 
theorem follows that K, defined in (|5.28p . has for each (xo,yo) sufficiently small a 
unique fixed point 

z = z(x ,y ) G E(R+,<5), 

and that the mapping [(xo,yo) z(xo,yo)] is continuous and vanishes at (0,0). 
By construction it follows that z = z(xo,yo) solves 

h Jt0J z = N(z QO + z + z) - N(Zoa), z(0) = (j}(z ) - 0(zoo). 

In summary, we have shown that for each Zq G SM.^ small enough, there exists 

(zoo, z, z) G Z x x E(M + , S) x E(R+, S) 

such that z = Zoo + z + z is the unique global solution of (|5.ip . In particular we 
have shown that for every zo G SAA~f small enough there exists a unique equilib- 
rium Zoo = ^oo(zo) such that the solution of (|5.1[) exists for alH > and converges 
to Zqo in SA4j at an exponential rate. 

(b) Now we consider the linearly unstable case and we first show that the equi- 
librium is unstable for the nonlinear equation (|5.1j) . Using the same notation as 
in part (a) we consider the system of equations 

h^ u z = N(z + z), z(0) = <f>(zo), 

z + LyS = ujz, z(0) = zq. 

Given a G R one verifies (by similar considerations as in [25l Proposition 10]) that 
there is a nondecreasing function r\ : R + — > R + such that r/(r) — > as r — > and 

\\e at N(z)l r{a) <r)(r)\\e at z\\ na) , e at zeE(a), (5.30) 
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whenever |^(t)|z < r for < t < a. Here a > is an arbitrary fixed number. For 
later use we note that 

E(a) M-L p ([0,a];X 7 ), (5.31) 
where the embedding constant is independent of a. 

Let a + be the collection of all positive eigenvalues of (— L 7 ). and let P + be 
the spectral projection related to the spectral set a + . Additionally, let P~ := 
I — P + and Xrf := P ± (X 1 ). Then X+ is finite dimensional and we obtain the 
decomposition 

X = X+®X~, L 1 = L+®L-. 
We note that a{— L+) = a + and <r(—L~) C [Rez < 0], where a(-L^) denotes 
the spectrum of (— Iy±), respectively. Let A* be the smallest positive eigenvalue 
of (— L+ ) and choose positive numbers k, [i such that [n — fi, n + [i] C (0, A*). We 
remind that the spectrum of (— L 7 ) consists of real eigenvalues, so that the strip 
[k — fi < Rez < k + fi] does not contain any spectral values of (— L 7 ). Therefore, 
there exists a constant M > 1 such that 

|e~ L ^| < Me*"-"' 1 , \e L ^\ < Me~ {K+tl)t , t>0. (5.32) 

Suppose now, by contradiction, that the equilibrium is stable for (|5.ip . Then 
for every r > there is a number S > such that (|5.1|) admits a global solution 
z e E(R + ) with |z(i)| < r for alH > whenever z e ^(0). 

In the following we will use the decomposition z — z + z, where (5, z) is the solution 
of the linear system (|5.29|) . (The function z = z + z is known, so that the first 
equation has a unique solution z. With z determined, z = z — z is the unique 
solution of the second equation.) The functions P ± z satisfy 

j t P ± z + L^P ± z = P ± luz 7 P ± z(0)=P ± z . (5.33) 

Next we shall show that P + z is given by the formula 

/oo 
e-^^-^P+uzdr, t>0. (5.34) 

Given any a > it follows from \P + z(t)\ x + < r that 

le- K ^ + zL p([ o,a ]; x+) < r {£ e ' Kptdt ) 1/P ^ °( K >Py ( 5 - 35 ) 
From the relation 

^- e - Kt P+z = f-K - L+)e- Kt P+z + e- Kt P + tuz, (5.36) 

(|535>(|536l) and (IOT1) follows 

\\e- Kt P+~z\\ x(a) < d(r + ||e- Kt z|| E(a) ), (5.37) 

with a universal constant C\. Here X(a) is defined as in (|5.16p . with the difference 
that R+ is replaced by the interval [0, a] and 6 — 0. We also recall that X+ is 
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finite dimensional, so that the spaces and D(L+) coincide (and therefore carry 
equivalent norms). From semigroup theory, maximal regularity, (|5.32p - (|5.33[) and 
(jOTI) follows 

\\e- Kt P-~z\\ x(a) < M(\P-z \ + \\e- Kt p-ujz\\ LpmahXy) ) 

<M(|p-^ |+C2||e- Kt f|| E(a) ). 

Combining (|5.37p - (|5.38l) results in 

|e- Kt % (a) < C 3 (r + \P-zo\ + |e- Kt z| K(a) ), (5.39) 

where C3 is a universal constant. Similarly as in part (a) we can infer from the 
equation for z that 

|e- Kt 5| E(Q) < c(\\e- Kt z\\ Ha) + ||e- Kt z|| E(Q) ), (5.40) 

and this implies 

||e- Kt z|| E(a) <C 4 {r + \P-2o\ + || e - Kt z|| E(a) ) (5.41) 

with C4 = c(l + C3). On the other hand we obtain from the equation for z and 
(1530)1 

|e- Kt z| E(a) < C(\4>(zo)\ + \\e- Kt N(S + z)\\ E(a) ) 

< C(\<P(z )\ +7 7 (r)(||e-' it z|| E(a) + ||e- Kt z|| E(a) )). 

If r is chosen small enough such that Crj(r) < 1/2 then 

\\e- Kt z\\ E(a} < 2C(\<j>(z )\ +v(r)\e- Kt z\ K{a) ). (5.42) 

We can, at last, combine (|5.41l) - (|5.42p to the result 

||e- K *z|| E(a) + ||e- K *z|| E(a) < C 5 (r + \P~zo\ + \4>{2o)\), (5.43) 

provided r is chosen small enough so that 2(l+C4)Ci~i(r) < 1/2. Since all estimates 
are independent of a we conclude that e~ Kt z G E(R + ). From (|5.43[) and Holder's 
inequality follows 

\e- L+ ^P+uz{T)\ x +dT 

a 00 1, , \ i/p' 
e w (t ~ T) dr) \e- Kr ojz\ Lv{VL+ .. x ^ < C\\ e - Kt z\\ m+) < 00, 

thus showing that the integral J t °° e~ L i P + ujz dr exists in for every t > 0. 
Moreover, its norm in grows no faster than the exponential function e Kt . It 
follows from the variation of parameters formula that 

/>oo />oo 

e L ^(P + z(t) + / e- L ^ {t - T) P + uJz{T) dr) = P + z + / e Lt < T P+ujz{t) dr, 
Jt Jo 

and the estimate 

/>oo 

\e L t t (P+z(t)+ e - L t {t - T) P + ujzdT)\ x+ < Me-^+^ir + Ce Kt ), t > 0, 
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then shows that P + zq + J °° e L i T P + ojz dr = 0. Thus the representation (|5.34[) 
holds as claimed. With this established, we obtain from Young's inequality for 
convolution integrals 

\\e- Kt P + m\\ Lp(R+ < ^M^le-^P+^I^^^). 

It then follows from (jOTj) and (j5~3BT) that 

\\e- Kt P + ~z\\x(m +) < C\\e- Kt z\\ m+) . (5.44) 

We can now imitate the estimates in (|5.38l) - (|5.42[) . with the interval [0, a] replaced 
by M. + , to conclude that 

\\e' Kt S\\E(R + ) + \\e- Kt z\\ m+) < C 6 (|P-fo| + \4>(zo)\). (5.45) 

This, in combination with (|5.32j) . (|5.34|) . and Holder's inequality, yields the esti- 
mate 

\P + zo\x+ < M / e-" T \e' KT P + ujz\ x + dr 

< C|| e - Kt P+ W z|| Lp(R+;X + ) < C(\P-~z \ + |0(z o )|). 

By decreasing S if necessary, we can assume that C|</>(5o)| < l/2(|-P + ^o| + l-P ^ol) 
for all 2 a e B s (0). (Recall that 0(0) = 0'(O) = 0.) Hence 

iP+Solz, <C 7 \P-z q \^, z eB s (0), (5.46) 

with a uniform constant CV, and this shows that cannot be stable for (|5 . 1 [) . 

It remains to show the last assertion of Theorem [521(b). For this we consider the 
projection P u = I — P c — P s which projects onto X™, the unstable subspace of 
Xj associated with the (finitely many) unstable eigenvalues. As in part (a) we 
will show that there exists an equilibrium Zqo such that any solution that stays in 
a small neighborhood of converges to = z co (zq) exponentially fast as t — » 00. 
Using the decomposition y = y s + y u , we obtain as in (a) the following system of 
equations: 

x = P c ljz, x(0)=x -Xoo, 
y s +£"y s = S' a (x 0O ,x,«), y s (0) = yg, (5.47) 
y u + L" y u = S u (xoo , x, z), y u (0) = y# , 
with 

Sj(Xao,X,z) = P J UJZ - 1p' j (x 00 + x)P c UJZ - Li f [tp j (x 00 + x) - Vj (*«>)]) 

where j 6 {s,u}, and where the functions ipj are defined similarly as in (|5 . 1 2[) . 
Suppose we have a global solution z € E(K + ) of ()5.1|1 with z(0) = z 6 <S.M 7 
which satisfies |z|j> < r, where r > is sufficiently small. By similar arguments 
as above (the presence of the function S u does not cause any principal difficulties) 
we infer that 

e- L ?(*- T >5 u (x OOJ x,z)dT ) t>0. (5.48) 
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For (xq,Yq,z) € Bi(ro,<5), with ro sufficiently small, we set 

/>oo 

Xoo := xo + / P c luz(t) dr, 
Jo 

/oo 
P c uz(t) dr, 

y. := f^ + £*,tr^ (S.(x OOJ x,«),yS) 
y„(i):=- / e-^-'^Xoo.x.zJdT. 



(5.49) 



As in part (a) we conclude that (|5.49|) admits for each (xo,yg,z) 6 Bi(ro,<5), with 
ro sufficiently small, a unique solution 

(xoo,x,y s ,y u ) = e(x ,yg,z) € X* x X c (M+,<5) x X s (R+,5) x X"(R+,<5). 

Following the arguments of part (a) then renders a solution 

3(x ,yo) = z 00 +x + i/j(x + x 00 ) - V(xoo) + y s +y« + z 

of (|5.ip with z = x + -0(x o ) + yg + yg + </>(x + V(*o) + yo + Yo)> where y ' is 
determined by 

/•OO 

yo = — / e L " T 5 u (x 001 x,g)dr. (5.50) 
Jo 

The solution 3(xo,yo) converges exponentially fast toward the equilibrium z^. In 
addition, we have shown that the initial value zq necessarily lies on the stable 
manifold belonging to z^, determined by the relation (|5.50[) . 

Due to uniqueness of (local) solutions to (|5.1[) . the solution 3(xo,yo) coincides 
with the given global solution z, and the proof of part (b) is now complete. □ 



Global existence and convergence. There are several obstructions against 
global existence for the Stefan problem 

• regularity: the norms of either u(t), T(t), and in addition [dc?„w(i)] m case 
7 = 0, become unbounded; 

• well-posedness: in case 7 = the well-posedness condition l(u) ^ may 
become violated; or u may become 0; 

• geometry: the topology of the interface changes; 
or the interface touches the boundary of fi; 

or the interface contracts to a point. 

Note that the compatibility conditions |/0( u )l + <J T~^ — in case 7 = 0, and 

(l(u) - Mu)j - *H)(Mu)] + aU) = y(u)ldd v u] 
in case 7 > are preserved by the semiflow. 
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Let (u, r) be a solution in the state manifold SM-y. By a uniform ball condition we 
mean the existence of a radius ro > such that for each t, at each point i £ T(t) 
there exist centers Xi € fli(t) such that B ro {xi) C fli and n B ro (xi) = {x\, 
i = 1,2. Note that this condition bounds the curvature of r(i), prevents it from 
shrinking to a point, from touching the outer boundary d£l, and from undergoing 
topological changes. 

With this property, combining the semiflow for (jl.ip with the Lyapunov func- 
tional and compactness we obtain the following result. 

Theorem 5.3. Let p > n + 2, a > 0, suppose ^,76 C 3 (Q, 00), d € C 2 (0, 00) such 
that either 7 = or j(u) > on (0, 00), and assume 

k(u) = — uip"(u) > 0, c?(u) > 0, we (0,oo). 

Suppose that [u, T) is a solution of in the state manifold SM^ on its maximal 
time interval [0,i*). Assume the following on [0,i*): 
W luWI^-v* + |rWl w 4-3/ P < M < 00; 

(ii) |[d9 l/ w(t)]| w/ 2-6/ P < M < 00 m case 7 = 0; 

(hi) > 1/M in case 7 ee 0; 

(iv) u(t) > 1/M; 

(v) r(i) satisfies a uniform ball condition. 

Then t* = 00, i.e. the solution exists globally. If its limit set contains a stable 
equilibrium (ttooiToo) € £, i.e. (^'(itoo) < 0, then it converges in SA4 7 to this 
equilibrium. On the contrary, if (u(t),T(t)) is a global solution in SA4 7 which 
converges to an equilibrium (u*,T*) with l(u*) =/= in case 7 = in SA4^ as 
t — > 00, then the properties (i)-(v) are valid. 

Proof. Assume that assertions (i)-(v) are valid. Then r([0,£*)) C Wp~ 3 ^ p (Q,r) is 
bounded, hence relatively compact in Wp 3 ^ p e (f2, r). (See Q3.7P for the definition 

ofw;(n,r).) 

Thus we may cover this set by finitely many balls with centers real analytic 
in such a way that dist^-3/p-e (r(t), S_y) < 5 for some j — j(t), t <E [0,i*). Let 
Jk = {t £ [0,t*) : j(t) = k}. Using for each k a Hanzawa-transformation S^, we 
see that the pull backs {u(t, -)oEk ■ t € Jk} are bounded in W^ _2y,p (f2\£fc), hence 
relatively compact in Wp - _e (fi\ £/-). Employing now Corollary 13. 101 we obtain 
solutions (m^T 1 ) with initial configurations (u(t),T(t)) in the state manifold on a 
common time interval, say (0,r], and by uniqueness we have 

(u\T),r 1 (r)) = (u(t + r),r(t + T)). 

Continuous dependence implies then relative compactness of (u(-),r(-)) in SM. 1 . 
In particular, = 00 and the orbit (w, T)(R + ) C SA4 7 is relatively compact. 
The negative total entropy is a strict Lyapunov functional, hence the limit set 
w(it, r) C SAi-y of a solution is contained in the set £ of equilibria. By compactness 
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u>(u, r) C SAi-y is non-empty, hence the solution comes close to £ , and stays there. 
Then we may apply the convergence result Theorem 15.21 The converse is proved 
by a compactness argument. □ 

Remark 5.4. We believe that the extra assumption (p'(uoo) < in Theorem 15.31 
can be replaced by (^'(itoo) 7^ 0. However, to prove this requires more technical 
efforts, and we refrain from doing this here. 
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